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SUMMARY 


This  document  is  the  final  report  for  an  investigation 
of  the  application  of  hierarchical  data  structures  to  geo¬ 
graphical  information  sys terns ,  under  Department  of  the  Army 
Contract  DAAK70-81-C-0059/P00007.  The  purposes  of  this 
investigation  were  twofold:  (1)  to  construct  a  geographic 
information  system  based  on  the  quadtree  hierarchical  data 
structure,  and  (2)  to  gather  statistics  to  allow  the  evalua¬ 
tion  of  the  usefulness  of  this  approach  to  geographic  infor¬ 
mation  system  organization.  To  accomplish  the  above  objec¬ 
tives,  a  database  was  built  that  contained  three  maps  sup¬ 
plied  under  the  terms  of  the  contract.  These  maps  described 
the  flood  plain,  elevation  contours,  and  landuse  classes  of 
a  region  in  California. 

This  study  report  presents  the  results  of  the  prelim¬ 
inary  investigation.  It  includes  analysis  of  the  merits  and 
deficiencies  of  the  various  approaches,  and  provides  recom¬ 
mendations  for  further  research. 
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l.\  Introduction 


~^This  project  is  concerned  with  the  applicability  of  a 
class  of  hierarchical  data  structures,  known  as  "quadtrees", 
to  the  representation  of  cartographic  data.  Section  2 
presents  a  tutorial  on  quadtree  data  structures.  Section  3 
describes  the  database  used,  and  the  process  of  digitizing 
and  editing  it.  Section  4  describes  the  process  of  quadtree 
encoding  of  the  data,  including  algorithms  and 
space/time/acreage  tables.  Section  5  discusses  region 
analysis  and  manipulations  using  quadtrees,  including  algor¬ 
ithms  and  tables  (time,  etc.).  The  algorithms  implemented 
include  set  theoretic  operations  on  regions,  point-in-region 
determination,  region  property  measurement,  and  construction 
of  submaps  and  merged  maps.  Section  6  presents  a  bibliogra¬ 
phy  on  quadtrees.  The  facilities  used  on  the  project  are 
described  in  the  Appendix. 

2.  Tutorial  on  quadtrees 

2.1.  Introduction 


In  gur  discussion  we  assume  that  a  region  is  a  subset 
of  a  2  by  2n  array  which  is  viewed  as  being  composed  of 
unit-square  pixels.  The  most  common  region  representations 
used  in  image  processing  are  the  binary  array  and  the  run 
length  representation  [1] .  The  binary  array  represents 
region  pixels  by  l's  and  non-region  Pixels  by  0's.  The  run 
length  representation  represents  each  row  of  the  binary 
array  as  a  sequence  of  runs  of  l's  alternating  with  runs  of 
0's. 


Boundaries  of  regions  are  often  specified  as  a  sequence 
of  unit  vectors  in  the  principal  directions.  This  represen¬ 
tation  is  termed  a  chain  code  [2J.  For  example,  letting  i 
represent  90  *  i  (i*0,l,2,3) ,  we  have  the  following  sequence 
as  the  chain  code  for  the  region  in  Figure  2.1a: 

0302352312330325160101030101 

Note  that  this  is  a  clockwise  code  which  starts  at  the  left¬ 
most  of  the  uppermost  border  points.  Chain  codes  yield  a 
compact  representation;  however,  they  are  somewhat  incon¬ 
venient  for  performing  operations  such  as  set  union  and 
intersection.  For  an  alternative  boundary  representation 
see  the  strip  trees  of  Ballard  [3] . 

Regions  can  also  be  represented  by  a  collection  of  max¬ 
imal  blocks  that  are  contained  in  the  given  region,  one 
such  trivial  representation  is  the  run  length  where  the 
blocks  are  1  by  m  rectangles.  A  more  general  representation 
treats  the  region  as  a  union  of  maximal  blocks  (of  l"s)  of  a 
given  shape.  The  medial  axis  transform  (MAT)  [4,5]  is  the 
set  of  points  serving  as  centers  of  these  blocks  and  their 
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Figure  2.1.  A  region,  its  maximal  blocks,  and  the  corresponding 

quadtree.  Blocks  in  the  region  are  shaded,  background 
blocks  are  blank.  Horizontal  lines  indicate  ropes. 
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corresponding  radii. 

The  quadtree  is  a  maximal  block  representation  in  which 
the  blocks  have  standard  sizes  and  positions  (i.e.,  powers 
of  two) .  It  is  an  approach  to  region  representation  which 
is  based  on  the  successive  subdivision  of  an  image  array 
into  quadrants.  If  the  array  does  not  consist  entirely  of 
l's  or  entirely  of  0's,  then  we  subdivide  it  into  quadrants, 
subquadrants,...  until  we  obtain  blocks  (possibly  single 
pixels)  that  consist  of  l's  or  of  0's,  i.e.,  they  are 
entirely  contained  in  the  region  or  entirely  disjoint  from 
it.  This  process  is  represented  by  a  tree  of  out  degree  4 
(i.e.,  each  non-leaf  node  has  four  sons)  in  which  the  root 
node  represents  the  entire  array.  The  four  sons  of  the  root 
node  represent  the  quadrants  (labeled  in  order  NW,  NE,  SW, 
SE) ,  and  the  leaf  nodes  correspond  to  those  blocks  of  the 
array  for  which  no  further  subdivision  is  necessary.  Leaf 
nodes  are  said  to  be  "black"  or  "white"  depending  on  whether 
their  corresponding  blocks  are  entirely  within  or  outside  of 
the  region  respectively.  All  non-leaf  nodes  are  said  to  be 
•gray".  Since  the  array  was  assumed  to  be  '  2n  by  2n,  the 
tree  height  is  at  most  n.  As  an  example.  Figure  2.1b  is  a 
block  decomposition  of  the  region  in  Figure  2.1a  while  Fig¬ 
ure  2.1c  is  the  corresponding  quadtree.  Each  quadtree  node 
is  implemented,  storage-wise,  as  a  record  with  six  fields. 
Five  fields  contain  pointers  to  the  four  sons  and  the  father 
of  a  node.  The  sixth  field  contains  type  information  such 
as  color,  etc.  Note  that  the  quadtree  representation  dis¬ 
cussed  here  should  not  be  confused  with  the  quadtree 
representation  of  two-dimensional  point  space  data  intro¬ 
duced  by  Finkel  and  Bentley  [6]  and  also  discussed  in  [7,8] 
and  improved  upon  in  [9] . 

The  quadtree  method  of  region  representation  is  based 
on  a  regular  decomposition.  it  has  been  employed  in  the 
domains  of  computer  graphics,  scene  analysis,  architectural 
design  [10] t  and  pattern  recognition.  In  particular, 
Warnock's  [10-13]  algorithm  for  hidden  surface  elimination 
is  based  on  such  a  principle — i.e.,  it  successively  subdi¬ 
vides  the  picture  into  smaller  and  smaller  squares  in  the 
process  of  searching  for  areas  to  be  displayed.  Application 
of  the  quadtree  to  image  representation  was  proposed  by 
Klinger  [14]  and  further  elaborated  upon  in  [15-20].  It  is 
relatively  compact  [15]  and  is  well  suited  to  operations 
such  as  union  and  intersection  [21-23],  and  detecting  vari¬ 
ous  region  properties  [15,21,22,  24].  Hunter's  Ph.D.  thesis 
[21,22,24],  in  the  domain  of  computer  graphics,  develops  a 
variety  of  algorithms  (including  linear  transformations)  for 
the  manipulation  of  a  quadtree  region  representation.  In 
[25-27]  variations  of  the  quadtree  are  applied  in  three 
dimensions  to  represent  solid  objects  and  in  [28]  to  more 
dimensions. 

There  has  been  much  work  recently  on  the 
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interchangeability  between  the  quadtree  and  other  tradi¬ 
tional  methods  of  region  representation.  Algorithms  have 
been  developed  for  converting  a  binary  array  to  a  quadtree 
[29] ,  run  lengths  to  a  quadtree  [30]  and  a  quadtree  to  run 
lengths  [31],  as  well  as  boundary  codes  to  a  quadtree  [32] 
and  a  quadtree  to  boundary  codes  [33] .  Work  has  also  been 
done  in  computing  geometric  properties  such  as  connected 
component  labeling  [34],  perimeter  [35],  Euler  number  [36], 
areas  and  moments  [23],  as  well  as  a  distance  transform 
[37,38].  In  addition,  the  quadtree  has  been  used  in  image 
processing  applications  such  as  shape  approximation  [39], 
edge  enhancement  [40] ,  image  segmentation  [41] ,  threshold 
selection  [42],  and  smoothing  [43]. 

2.2.  Preliminaries 


In  the  quadtree  representation,  by  virtue  of  its  tree¬ 
like  nature,  most  operations  are  carried  out  by  techniques 
which  traverse  the  tree.  In  fact,  many  of  the  operations 
that  we  describe  can  be  characterized  as  having  two  basic 
steps.  The  first  step  either  traverses  the  quadtree  in  a 
specified  order  or  constructs  a  quadtree.  The  second  step 
performs  a  computation  at  each  node  which  often  makes  use  of 
its  neighboring  nodes,  i.e.,  nodes  representing  image  blocks 
that  are  adjacent  to  the  given  ncde's  block.  For  examples, 
see  [30-38].  Frequently,  these  two  steps  are  performed  in 
parallel . 

In  general,  it  is  preferable  to  avoid  having  to  use 
position  (i.e.,  coordinates)  and  size  information  when  mak¬ 
ing  relative  transitions  (i.e.,  locating  neighboring  nodes) 
in  the  quadtree  since  they  involve  computation  (rather  than 
simply  chasing  links)  and  are  clumsy  when  adjacent  blocks 
are  of  different  sizes  (e.g.,  when  a  neighboring  block  is 
larger) .  similarly,  we  do  not  assume  that  there  are  links 
from  a  node  to  its  neighbors,  because  we  do  not  want  to  use 
links  in  excess  of  four  links  from  a  non-leaf  node  to  its 
sons  and  the  link  from  a  non-root  node  to  its  father.  Such 
techniques,  described  in  [44],  are  used  in  [30-38]  and 
result  in  algorithms  that  only  make  use  of  the  existing 
structure  of  the  tree.  This  is  in  contrast  with  the  methods 
of  Klinger  and  Rhodes  [19]  which  make  use  of  size  and  posi¬ 
tion  information,  and  those  of  Hunter  and  Steiglitz  [21, 
22,24]  which  locate  neighobrs  through  the  use  of  explicit 
links  (termed  nets  and  ropes) . 

Locating  neighbors  in  a  given  direction  is  quite 
straightforward.  Given  a  node  corresponding  to  a  specific 
block  in  the  image,  its  neighbor  in  a  particular  direction 
(horizontal  or  vertical)  is  determined  by  locating  a  common 
ancestor.  For  example,  if  we  want  to  find  an  eastern  neigh¬ 
bor,  the  common  ancestor  is  the  first  ancestor  node  which  is 
reached  via  its  NW  or  SW  son.  Next,  we  retrace  the  path 
from  the  common  ancestor,  but  making  mirror  image  moves 
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about  the  appropriate  axis,  e.g.,  to  £ind  an  eastern  or 
western  neighbor,  the  mirror  images  of  NE  and  SE  are  NW  and 
SW,  respectively.  For  example,  the  eastern  neighbor  of  node 
32  in  Figure  2.1c  is  node  33.  It  is  located  by  ascending 
the  tree  until  the  common  ancestor,  H,  is  found.  This 
requires  going  through  a  SE  link  to  reach  L  and  a  NW  link  to 
reach  H.  Node  33  is  now  reached  by  backtracking  along  the 
previous  path  with  the  appropriate  mirror  image  moves  (i.e., 
going  through  a  NE  link  to  reach  M  and  a  SW  link  to  reach 
33)  . 


In  general,  adjacent  neighbors  need  not  be  of  the  same 
size.  If  they  are  larger,  then  only  a  part  of  the  path  to 
the  common  ancestor  is  retraced.  If  they  are  smaller,  then 
the  retraced  path  ends  at  a  "gray"  node  of  equal  size.  Thus 
a  "neighbor"  is  correctly  defined  as  the  smallest  adjacent 
leaf  whose  corresponding  block  is  of  greater  than  or  equal 
size.  If  no  such  node  exists,  then  a  gray  node  of  equal 
size  is  returned.  Note  that  similar  techniques  can  be  used 
to  locate  diagonal  neighbors  (i.e.,  nodes  corresponding  to 
blocks  that  touch  the  given  node's  block  at  a  corner).  For 
example,  node  20  in  Figure  2.1c  is  the  NW  neighbor  of  node 
22.  For  more  details,  see  [44]. 

In  contrast  with  our  neighbor  finding  methods  is  the 
use  of  explicit  links  from  a  node  to  its  adjacent  neighbors 
in  the  horizontal  and  vertical  directions  reported  in 
[21,22,24] .  This  is  achieved  through  the  use  of  adjacency 
trees,  "ropes,"  and  "nets."  An  adjacency  tree  exists  when¬ 
ever  a  leaf  node,  say  X,  has  a  GRAY  neighbor,  say  Y,  of 
equal  size.  In  such  a  case,  the  adjacency  tree  of  X  is  a 
binary  tree  rooted  at  Y  whose  nodes  consist  of  all  sons  of  Y 
(BLACK,  WHITE,  and  GRAY)  that  are  adjacent  to  X.  For  exam¬ 
ple,  for  node  16  in  Figure  2.1,  the  western  neighbor  is  GRAY 
node  F  with  an  adjacency  tree  as  shown  in  Figure  2.2.  A  rope 
is  a  link  between  adjacent  nodes  of  equal  size  at  least  one 
of  which  is  a  leaf  node.  For  example,  in  Figure  2.1,  there 
exists  a  rope  between  node  16  and  nodes  G,  17,  H,  and  F. 
Similiarly,  there  exists  a  rope  between  node  37  and  nodes  N 
and  N;  however,  there  does  not  exist  a  rope  between  node  L 
and  nodes  M  and  N. 

The  algorithm  for  finding  a  neighbor  using  a  roped 
quadtree  is  quite  simple.  We  want  a  neighbor,  say  Y,  on  a 
given  side,  say  D,  of  a  block,  say  X.  If  there  is  a  rope 
from  X  on  side  D,  then  it  leads  to  the  desired  neighbor.  If 
no  such  rope  exists,  then  the  desired  neighbor  must  be 
larger.  In  such  a  case,  we  ascend  the  tree  until  encounter¬ 
ing  a  node  having  a  rope  on  side  D,  that  leads  to  the 
desired  neighbor.  In  effect,  we  have  ascended  the  adjacency 
tree  of  Y.  For  example,  to  find  the  eastern  neighbor  of  node 
21  in  Figure  2.1,  we  ascend  through  node  J  to  node  F,  which 
has  a  rope  along  its  eastern  side  leading  to  node  16. 


lqure  2.2. 


rtdjacency  tree  for  the  western  neighbor  oi  none  it> 
in  Figure  2  .i  . 


2  .<*  . 


uiocKg  i.  ana  w  ending  at  a  common  corner  . 
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At  times  it  is  not  convenient  to  ascend  nodes  searching 
for  ropes.  A  data  structure  named  a  net  is  used  [21,  22, 
24]  to  obviate  this  step  by  linking  all  leaf  nodes  to  their 
neighbors  regardless  of  their  size.  Thus  in  the  previous 
j  example  there  would  be  a  direct  link  between  nodes  21  and  16 

along  the  eastern  side  of  node  21.  The  advantage  of  ropes 
and  nets  is  that  the  number  of  links  that  must  be  traversed 
is  reduced.  However,  the  disadvantage  is  that  the  storage 
requirements  are  considerably  increased  since  many  addi¬ 
tional  links  are  necessary.  In  contrast,  our  methods  are 
I  implemented  by  algorithms  that  make  use  of  the  existing 

structure  of  the  tree  —  i.e.,  four  links  from  a  nonleaf 
node  to  its  sons,  and  a  link  from  a  nonroot  node  to  its 
father. 


2.3.  Conversion 

2.3.1.  Quadtrees  and  Arrays 

The  definition  of  a  quadtree  leads  naturally  to  a  "top 
down”  quadtree  construction  process.  This  may  lead  to 
excessive  computation  because  the  process  of  examining 
whether  a  quadtrant  contains  all  l's  or  all  0's  may  cause 
certain  parts  of  the  region  to  be  examined  repeatedly  by 
virtue  of  being  composed  of  a  mixture  of  l's  and  O's. 
Alternatively,  a  ” bottom-up”  method  may  be  employed  which 
scans  the  picture  in  the  sequence 


1 

2 

5 

6 

17 

18 

21 

22 

3 

4 

7 

8 

19 

20 

23 

24 

9 

10 

13 

14 

25 

26 

29 

30 

11 

12 

15 

16 

27 

28 

31 

32 

33  ... 

where  the  numbers  indicate  the  sequence  in  which  the  pixels 
are  examined.  As  maximal  blocks  of  0's  or  l*s  are 
discovered,  corresponding  leaf  nodes  are  added  along  with 
the  necessary  ancestor  nodes.  This  is  done  in  such  a  way 
that  leaf  nodes  are  never  created  until  they  are  known  to  be 
maximal.  Thus  there  is  never  a  need  to  merge  four  leaves  of 
the  same  color  and  change  the  color  of  their  common  parent 
from  gray  to  white  or  black  as  is  appropriate.  See  [29]  for 
the  details  of  such  an  algorithm  whose  execution  time  is 
proportional  to  the  number  of  pixels  in  the  image. 

If  it  is  necessary  to  scan  the  picture  row  by  row 
(e.g.,  when  the  input  is  a  run  length  coding)  the  quadtree 
construction  process  is  somewhat  more  complex,  we  scan  the 
picture  a  row  at  a  time.  For  odd-numbered  rows,  nodes 
corresponding  to  the  pixel  or  run  values  are  added  for  the 
pixels  and  attempts  are  made  to  discover  maximal  blocks  of 
0*8  or  l's  whose  size  depends  on  the  row  number  (e.g.,.  when 
processing  the  fourth  row,  maximal  blocks  of  maximum  size 
4- by- 4  can  be  discovered) .  In  such  a  case  merging  is  said 
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to  take  place.  See  [30]  for  the  details  of  an  algorithm 
that  constructs  a  quadtree  from  a  row  by  row  scan  such  that 
at  any  instance  of  time  a  valid  quadtree  exists.  This 
algorithm  has  an  execution  time  that  is  proportional  to  the 
number  of  pixels  in  the  image. 

Similarly,  for  a  given  quadtree  we  can  output  the 
corresponding  binary  picture  by  traversing  the  tree  in  such 
a  way  that  for  each  row  the  appropriate  blocks  are  visited 
and  a  row  of  0's  or  l's  is  output.  Zn  essence,  we  visit 
each  quadtree  node  once  for  each  row  that  intersects  it 
(i.e.,  a  node  corresponding  to  a  block  of  size  2 K  by  2K  is 
visited  2K  times).  For  the  details  see  [31]  where  an  algor¬ 
ithm  is  described  whose  execution  time  depends  only  on  the 
number  of  blocks  of  each  size  that  comprise  the  image  -  not 
on  their  paticular  configuration. 

2.3.2.  Quadtrees  and  borders 

In  order  to  determine,  for  a  given  leaf  node  M  of  a 
quadtree,  whether  the  corresponding  block  is  on  the  border, 
we  must  visit  the  leaf  nodes  that  correspond  to  4-adjacent 
blocks  and  check  whether  they  are  black  or  white.  For  exam¬ 
ple,  to  find  M's  right  hand  neighbor  in  Figure  2.3,  we  use 
the  neighbor  finding  techniques  outlined  in  Section  2.2.  If 
the  neighbor  is  a  leaf  node,  then  its  block  is  at  least  as 
larga  as  that  of  M  and  so  it  is  M's  sole  neighbor  to  the 
right.  Otherwise,  the  neighbor  is  the  root  of  a  subtree 
whose  leftmost  leaf  nodes  correspond  to  M's  right-hand 
neighbors.  These  nodes  are  found  by  traversing  that  sub¬ 
tree. 


Let  M,N  in  Figure  2.3  be  black  and  white  leaf  nodes 
whose  associated  blocks  are  4-adjacent.  Thus  the  pair  M,N 
defines  a  common  border  segment  of  length  2K  (2K  is  the 
minimum  of  the  side  lengths  of  M  arid  N)  which  ends  at  a 
corner  of  the  smaller  of  the  two  blocks  (they  may  both  end 
at  a  common  point  as  in  Figure  2.4).  In  order  to  produce  a 
boundary  code  representation  for  a  region  in  the  image  we 
must  determine  the  next  segment  along  the  border  whose  pre¬ 
vious  segment  lay  between  M  and  N.  This  is  achieved  by 
locating  the  other  leaf  P  whose  block  touches  the  end  of  the 
segment  between  M  and  N.  If  the  M,N  segment  ends  at  a 
corner  of  both  M  and  N,  then  we  must  find  the  other  leaf  R 
or  leaves  P,Q  whose  blocks  touch  that  corner  (see  Figure 
2.4)  Again,  this  can  be  accomplished  by  using  neighbor  find¬ 
ing  techniques  as  outlined  in  Section  2.2. 

For  the  non-common  corner  case,  the  next  border  segment 
is  the  common  border  defined  by  M  and  P  if  P  is  white,  or 
the  common  border  defined  by  N  and  p  if  p  is  black,  in  the 
common  corner  case,  the  pair  of  blocks  defining  the  next 
border  segment  is  determined  exactly  as  in  the  standard 
"crack  following"  algorithm  [45]  for  traversing  region 
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borders.  This  process  is  repeated  until  we  re-encounter  the 
block  pair  M,N.  At  this  point  the  entire  border  has  been 
traversed.  The  successive  border  segments  constitute  a  4- 
direction  chain  code,  broken  up  into  segments  whose  lengths 
are  sums  of  powers  of  two.  The  time  required  for  this  pro¬ 
cess  is  on  the  order  of  the  number  of  border  nodes  times  the 
tree  height.  For  more  details  see  [33]  . 

Using  the  methods  described  in  the  last  two  paragraphs, 
we  can  traverse  the  quadtree,  find  all  borders,  and  generate 
their  codes.  During  this  process,  we  mark  each  border  as  we 
follow  it,  so  that  it  will  not  be  followed  again  from  a  dif¬ 
ferent  starting  point.  Note  that  the  marking  process  is 
complicated  by  the  fact  that  a  node's  block  may  be  on  many 
different  borders. 

In  order  to  generate  a  quadtree  from  a  set  of  4- 
direction  chain  codes  we  use  a  two-step  process.  First,  we 
trace  the  boundary  in  a  clockwise  direction  and  construct  a 
quadtree  whose  black  leaf  nodes  are  of  a  size  equal  to  the 
unit  code  length.  All  the  black  nodes  correspond  to  blocks 
on  the  interior  side  of  the  boundary.  All  remaining  nodes 
are  left  uncolored.  Second,  all  uncolored  nodes  are  set  to 
black  or  white  as  appropriate.  This  is  achieved  by  travers¬ 
ing  the  tree,  and  for  each  uncolored  leaf  node,  examining 
its  neighbors.  The  node  is  colored  black  unless  any  of  its 
neighbors  is  white  or  is  black  with  a  border  along  the 
shared  boundary.  At  any  stage,  merging  occurs  if  the  four 
rows  of  a  non-leaf  node  are  leaves  having  the  same  color. 
The  details  of  the  algorithm  are  given  in  [32].  The  time 
required  is  proportional  to  the  product  of  the  perimeter 
(i.e.,  the  4-direction  chain  code  length)  and  the  tree 
height. 


2.3.31.  Quadtrees  of  derived  sets 

Let  S  be  the  set  of  l's  in  a  given  binary  array,  and 
let  S  be  the  complement  of  S.  The  quadtree  of  the  comple¬ 
ment  of  S  is  the  same  as  that  of  S,  with  black  leaf  nodes 
changed  to  white  and  vice  versa.  To -.get  the  quadtree  of  the 
union  of  S  and  T  from  those  of  S  and  T,  we  traverse  the  two 
trees  simultaneously.  Where  they  agree,  the  new  tree  is  the 
same  and  if  the  two  nodes  are  gray,  then  their  subtrees  are 
traversed.  If  S  has  a  gray  (*nonleaf)  node  where  T  has  a 
black  node,  the  new  tree  gets  a  black  node;  if  T  has  a  white 
node  there,  we  copy  the  subtree  of  s  at  that  gray  node  into 
the  new  tree.  If  S  has  a  white  node,  we  copy  the  subtree  of 
T  at  the  corresponding  node.  The  algorithm  for  the  inter¬ 
section  of  S  and  T  is  exactly  analogous,  with  the  roles  of 
black  and  white  reversed.  The  time  required  for  these 
algorithms  is  proportional  to  the  number  of  nodes  in  the 
smaller  of  the  two  trees  [23]. 


10 


2.3.4.  Skeletons  and  medial  axis  transforms 

The  medial  axis  of  a  region  is  a  subset  of  its  points 
each  of  which  has  a  distance  from  the  complement  of  the 
region  (using  a  suitably  defined  distance  metric)  which  is  a 
local  maximum.  The  medial  axis  transform  (MAT)  consists  of 
the  set  of  medial  axis  or  "skeleton"  points  and  their  asso¬ 
ciated  distance  values.  The  quadtree  representation  may  be 
rendered  even  more  compact  by  the  use  of  a  skeleton-like 
representation.  Recall  that  a  quadtree  is  a  set  of  disjoint 
maximal  square  blocks  having  sides  whose  lengths  are  powers 
of  2.  We  define  a  quadtree  skeleton  to  be  a  set  of  maximal 
square  blocks  having  sides  whose  lengths  are  sums  of  powers 
of  two.  The  maximum  value  (i.e.,  "chessboard")  distance 
metric  [45]  is  the  most  appropriate  for  an  image  represented 
by  a  quadtree.  See  [37]  for  the  details  of  its  computation 
for  a  quadtree;  see  also  [38]  for  a  different  quadtree  dis¬ 
tance  transform.  A  quadtree  medial  axis  transform  (QMAT)  is 
a  quadtree  whose  black  nodes  correspond  to  members  of  the 
quadtree  skeleton  while  all  remaining  leaf  nodes  are  white. 
The  QMAT  has  several  important  properties.  First,  it 
results  in  a  partition  of  the  image  into  a  set  of  possibly 
non-disjoint  squares  having  sides  whose  lengths  are  sums  of 
powers  of  two  rather  than,  as  is  the  case  with  quadtrees,  a 
set  of  disjoint  squares  having  sides  of  lengths  which  are 
powers  of  two.  Second,  the  QMAT  is  more  compact  than  the 
quadtree  and  has  a  decreased  shift  sensitivity.  See  [46] 
for  the  details  of  a  quadtree  to  QMAT  conversion  algorithm 
whose  execution  time  is  on  the  order  of  the  number  of  nodes 
in  the  tree. 

i  2.4.  Property  measurement 
2.4.1,.  Connected  component  labeling 

Traditionally,  connected  component  labeling  is  achieved 
by  scanning  a  binary  array  row  by  row  from  left  to  right  and 
labeling  adjacencies  that  are  discovered  to  the  right  and 
downward.  During  this  process  equivalences  will  be  gen¬ 
erated.  A  subsequent  pass  merges  these  equivalences  and 
updates  the  labels  of  the  affected  pixels.  In  the  case  of 
the  quadtree  representation  we  also  scan  the  image  in  a 
sequential  manner.  However,  the  sequence's  order  is  dic¬ 
tated  by  the  tree  structure  -  i.e.,  we  traverse  the  tree  in 
postorder.  Whenever  a  black  leaf  node  is  encountered  all 
black  nodes  that  are  adjacent  to  its  south  and  east  sides 
are  also  visited  and  are  labeled  accordingly.  Again, 
equivalences  generated  during  this  traversal  are  subse¬ 
quently  merged  and  a  tree  traversal  is  used  to  update  the 
labels.  The  interesting  result  is  that  the  algorithm's  exe¬ 
cution  time  is  proportional  to  the  number  of  pixels.  An 
analgous  result  is  described  in  the  next  section.  See  [34] 
for  the  details  of  an  algorithm  that  labels  connected  com¬ 
ponents  in  time  on  the  order  of  the  number  of  nodes  in  the 
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tree  plus  the  product  of  B'log  B  where  B  is  the  number  of 
black  leaf  nodes. 

2.4.2.  Component  counting  and  genus  computation 


Once  the  connected  components  have  been  labeled,  it  is 
trivial  to  count  them,  since  their  number  is  the  same  as  the 
number  of  inequivalent  labels.  We  will  next  describe  a 
method  of  determining  the  number  of  components  minus  the 
number  of  holes  by  counting  certain  types  of  local  patterns 
in  the  array;  this  number,  g,  is  known  as  the  genus  or  Euler 
number  of  the  array. 

Let  V  be  the  number  of  l's,  E  the  number  of  horizon¬ 
tally  adjacent  pairs  of  l's  (i.e.,  11)  and  vertically  adja¬ 
cent  pairs  of  l's,  and  F  the  number  of  two  by  two  arrays  of 
l's  in  the  array;  it  is  well  known  [45)  that  g*v-E+F.  This 
result  can  be  generalized  to  the  case  where  the  array  is 
represented  by  a  quadtree  [36].  In  fact,  let  V  be  the 
number  of  black  leaf  nodes;  E  the  number  of  pairs  of  such 
nodes  whose  blocks  are  horizontally  or  vertically  adjacent; 
and  F  the  number  of  triples  or  quadruples  of  such  nodes 
whose  blocks  meet  at  and  surround  a  common  point  (see  Figure 
2.5) .  Then  g«V-E+F.  These  adjacencies  can  be  found  (see 
section  2.3.2)  by  traversing  the  tree;  the  time  required  is 
on  the  order  of  the  number  of  nodes  in  the  tree. 

2.4.3.  Area  and  moments 

The  area  of  a  region  represented  by  a  quadtree  can  be 
obtained  by  summing  the  areas  of  the  black  leaf  nodes,  i.e.* 
counting  4n  for  each  such  node  that  represents  a  2n  by  2n 
block.  Similarly,  the  first  x  and  y  moments  of  the  region 
relative  to  a  given  origin  can  be  computed  by  summing  the 
first  moments  of  these  blocks;  note  that  we  know  the  posi¬ 
tion  (and  size)  of  each  block  from  the  coordinates  of  its 
leaf  in  the  tree.  Knowing  the  area  and  the  first  moments 
gives  us  the  coordinates  of  the  centroid,  and  we  can  then 
compute  central  moments  relative  to  the  centroid  as  the  ori¬ 
gin.  The  time  required  for  any  of  these  computations  is 
proportional  to  the  number  of  nodes  in  the  tree.  Further 
details  on  moment  computation  from  quadtrees  can  be  found  in 
123). 

2.4.4.  perimeter 

An  obvious  way  of  obtaining  the  perimeter  of  a  region 
represented  by  a  quadtree  is  to  simply  traverse  its  border 
and  sum  the  number  of  steps.  However,  there  is  no  need  to 
traverse  the  border  segments  in  order.  Instead,  we  use  a 
method  which  traverses  the  tree  in  postorder  and  for  each 
black  leaf  node  examines  the  colors  of  its  neighbors  on  its 
four  sides.  For  each  white  neighbor  the  length  of  the 
corresponding  border  segment  is  included  in  the  perimeter. 


See  [35]  for  the  details  of  such  an  algorithm  which  has  exe¬ 
cution  time  proportional  to  the  number  of  nodes  in  the  tree. 
An  even  better  formulation  is  reported  in  [47]  which  gen¬ 
eralizes  the  concept  of  perimeter  to  n  dimensions. 

2S.5.  Concluding  remarks 

We  have  briefly  sketched  algorithms  for  accomplishing 
traditional  region  processing  operations  by  use  of  the  quad¬ 
tree  representation.  Many  of  the  methods  used  on  the  pixel 
level  carry  over  to  the  quadtree  domain  (e.g.,  connected 
component  labeling,  genus,  etc.).  Because  of  its  compact¬ 
ness,  the  quadtree  permits  faster  execution  of  these  opera¬ 
tions.  Often  the  quadtree  algorithms  require  time  propor¬ 
tional  to  the  number  of  blocks  in  the  image,  independent  of 
their  size. 


The  quadtree  data  structure  requires  storage  for  the 
various  links.  However,  use  of  neighbor  finding  techniques 
rather  than  ropes  a  la  Hunter  [21,  22,  24]  is  a  compromise. 
In  fact,  experimental  results  discussed  in  the  data  analysis 
segment  of  this  report  show  that  the  extra  storage  cost  of 
ropes  is  not  justified  by  the  resulting  minor  decrease  in 
execution  time.  This  is  because  the  average  number  of  links 
traversed  by  neighbor  finding  methods  is  3.5  in  contrast 
with  1.5  for  ropes.  Nevertheless,  there  is  a  possibility 
that  the  quadtree  may  not  be  efficient  spacewise.  For  exam¬ 
ple,  a  checkerboard-like  region  does  not  lead  to  economy  of 
space.  The  space  efficiency  of  the  quadtree  is  analyzed  in 
[48].  Some  savings  can  be  obtained  by  normalizing  the  quad¬ 
tree  [49,50]  as  is  also  possible  by  constructing  a  forest  of 
quadtrees  [51]  to  avoid  large  regions  of  WHITE.  Storage  can 
also  be  saved  by  using  a  locational  code  for  all  BLACK 
blocks  [52].  Gray  level  quadtrees  using  a  sequence  of  array 
codes  to  economize  on  storage  are  reported  in  [53] . 

The  quadtree  is  especially  useful  for  point  in  polygon 
operations  as  well  as  for  query  operations  involving  image 
overlays  and  set  operations.  The  hierarchical  nature  enables 
one  to  use  image  approximations,  in  particular,  a  breadth- 
first  transmission  of  an  image  yields  a  successively  finer 
image  yet  enabling  the  user  to  have  a  partial  image.  Thus 
the  quadtree  could  be  used  in  browsing  through  a  large  image 
database. 

Quadtrees  constitute  an  interesting  alternative  to  the 
standard  methods  of  digitally  representing  regions.  Their 
chief  disadvantage  is  that  they  are  not  shift-invariant;  two 
regions  differing  only  by  a  translation  may  have  quite  dif¬ 
ferent  quadtrees  (but  see  [46]).  Thus  shape  matching  from 
quadtrees  is  not  straightforward.  Nevertheless,  in  other 
respects,  they  have  many  potential  advantages.  They  provide 
a  compact  and  easily  constructed  representation  from  which 
standard  region  properties  can  be  efficiently  computed.  In 
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effect,  they  are  "variable-resolution  arrays"  in  which 
detail  is  represented  only  when  it  is  available,  without 
requiring  excessive  storage  for  parts  of  the  image  where 
detail  is  missing.  Their  variable-resolution  property  is 
superior  to  trees  based  on  a  hexagonal  decomposition  [54] 
in  that  a  square  can  be  repeatedly  decomposed  into  smaller 
squares  (as  can  be  done  for  triangles  as  well  [55])  whereas 
once  the  smallest  hexagon  has  been  chosen  it  can  not  be 
further  decomposed  into  smaller  hexagons.  Note  that  the 
variance  of  resolution  only  applies  to  the  area.  For  an 
application  of  the  quadtree  concept  to  borders,  as  well  as 
area,  see  the  line  quadtree  of  [56]. 


15 


3.  Database,  digitization,  and  editing 
3.1.  Procedures  and  results 

The  data  supplied  by  ETL  consisted  of  three  map  over¬ 
lays  (Pigures  3.1-3)  representing  land  use  classes,  terrain 
elevation  contours,  and  flood  plain  boundaries  for  a  small 
area  of  Northern  California.  These  overlays  are  shown,  at  a 
reduced  scale,  in  the  figures  attached  to  this  section.  In 
the  case  of  the  elevation  contours,  only  those  at  multiples 
of  100  feet  were  to  be  digitized,  and  for  all  three  over¬ 
lays,  only  the  portions  bounded  by  the  fiducial  marks. 

Conversion  of  the  data  to  machine-readable  form  was 
carried  out  as  follows:  Each  overlay  was  superimposed  on  a 
grid  (graph  paper,  20  boxes  to  the  inch).  The  boundaries  to 
be  digitized  were  followed  by  hand  and  marked  on  a  second 
sheet  of  graph  paper.  Every  box  on  the  original  graph  was 
copied  onto  a  2-by-2  block  of  boxes  on  the  second  sheet. 
This  yielded  increased  resolution  and  also  separated  boun¬ 
dary  lines  which  on  the  original  graph  would  have  been  in 
adjacent  boxes.  This  graph  was  then  hand  chain-coded  and 
the  chain-codes  were  typed  into  the  computer  (see  the 
description  of  the  program  "mkbin"  for  a  definition  of  the 
chain-code  used) . 

A  binary  array  was  created  for  each  of  the  three 
overlays  in  which  the  pixels  that  were  on  a  boundary  in  the 
original  overlay  are  represented  by  a  value  of  1,  and  all 
other  pixels  are  represented  by  a  value  of  0.  A  connected 
component  labeling  program  was  then  applied  to  this  array 
yielding  an  array  in  which  the  pixels  in  each  connected 
region  have  a  unique  label.  (Pixels  of  value  1  were 
regarded  as  connected  even  if  they  were  only  diagonally 
adjacent,  whereas  pixels  of  value  0  were  regarded  as  con¬ 
nected  only  if  they  were  horizontally  or  vertically  adja¬ 
cent.)  A  lookup  table  was  then  created  to  convert  these 
labels  to  a  consistent  label  set  in  which  all  regions  of  a 
given  land-use  class, or  all  regions  between  a  given  pair  of 
elevation  contours,  had  the  same  label.  At  this  time,  all 
polygons  on  the  landuse  map  which  either  had  no  label  or  for 
which  the  label  was  unclear  were  placed  in  a  special  landuse 
class  "unk". 

The  final  data  preparation  task  was  to  remove  the 
boundary  lines  separating  the  regions  (those  pixels  given  a 
value  of  1  in  the  binary  array) .  This  was  uniformly  done  by 
assigning  to  each  boundary  pixel  the  label  of  its  right-hand 
neighbor,  or  if  this  was  also  a  boundary  pixel,  the  label  of 
its  neighbor  in  the  row  above.  The  three  digital  maps  (one 
per  overlay)  resulting  from  this  processing  were  450  pixels 
high  by  400  wide,  partitioned  into  labelled  regions  with  no 
"black"  boundary  lines  separating  them. 
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The  resulting  maps  were  then  quadtree  encoded  using 
the  quadtree  building  algorithm  described  in  the  Section  4. 
For  each  map,  a  multi-color  quadtree  was  built  giving  every 
region  in  the  map  a  unique  label.  A  multi-color  quadtree 
refers  to  a  quadtree  in  which  the  leaf  nodes  can  have  dif¬ 
ferent  colors.  Thus  a  multi-color  quadtree  is  an  extension 
of  the  black  and  white  (or  binary)  quadtree  that  is  dis¬ 
cussed  in  section  2.  Certain  operations,  for  example  union 
and  intersection,  are  not  defined  in  terms  of  multi-color 
quadtrees,  but  rather  are  defined  in  terms  of  binary  quad¬ 
trees  that  are  derived  from  the  multi-color  quadtrees  by 
considering  one  of  the  colors  as  black  (this  is  usually  the 
color  of  the  object  of  interest)  and  the  other  colors  to  be 
white.  In  addition  to  the  above  three  multi-color  quad¬ 
trees,  a  quadtree  was  built  for  each  land-use  class  or 
elevation  level  (i.e.  regions  in  the  class  or  elevation  are 
labeled,  all  other  regions  are  white) .  Programs  were  then 
written  to  manipulate  and  display  quadtrees,  as  well  as  cal¬ 
culate  region  properties  and  compute  set  theoretic  opera¬ 
tions  on  the  trees.  These  programs  are  described  in  later 
sections. 

The  hand  digitization  process  took  approximately  100 
manhours,  including  both  planning  and  implementation.  This 
time  could  be  greatly  shortened  by  using  coordinate  digitiz¬ 
ing  equipment.  Editing  of  the  hand-input  data  was  carried 
out  by  visual  inspection  of  the  resulting  regions  to  verify 
that  there  were  no  gaps  or  overlaps.  This  process,  together 
with  a  few  hand  corrections  of  touching  lines,  took  at  most 
20  manhours. 

Figures  3.4-3.38  show  the  components  of  each  land  use 
class.  Figures  3.39-3.49  show  the  components  of  each  eleva¬ 
tion  level,  and  Figure  3.50  shows  the  three  components  of 
the  flood  plain  map. 


fiaure  3. 3.  Yne  13  components  or  tne  lana-use  class  ACP . 


fiaure  3.o.  'i'ne  b  components  or  cue  land-use  class  ah . 


figure  j.7.  Tne  i  component  of  tne  land-use  class  HKt . 


J.io.  me  2  components  or  tne  lana-use  class  t . 


figure  3.14.  i’ne  s  components  o£  tne  lanci-us-^  class  ru. 
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riaura  3.17.  The  2  components  or  tne  land-use  class  ucu 


r'iaure  3.20.  The  2  components  oi  tne  lana-use  class  ucm. 


r'igure  3.21  Tne  2  components  ox  tne  lana-use  class  Ux.£> . 
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figure  The  2  components  or  tne  land-use  cl 
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figure  3.2o.  The  3  components  of  tne  lana-use  class  uoo . 


igure  3.2y.  The  2  components  of  the  lana-use  class  UuP . 


Figure  3.30.  Tne  2  components  of  tne  lana-use  class  Uuv . 


fcigure  3.32.  The  2h  components  ot  the  lanci-use  ciass  u«i» . 


Figure  3.33.  Tne  2  components  ot  the  land-use  class  UUi>. 
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r'zgure  3.3©.  'i’ne  2  components  ot  the  lana-use 


figure  3.37.  The  i  component  of  the  land-use  class  . 
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nqure  3.42.  Vue  1j  components  o t  tne  4tn  elevation  level. 

(300  -  400  ft.) 


Ei 


figure  3.4i.  Tne  d  components  ot  tne  ?tn  elevation  level. 

(600  -  700  ft.) 
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rigure  3.47.  Tne  o  components  of  tike  yth  elevation  level. 

(800  -  900  ft.) 


components  of  tne  ilood-plain  i^ap. 


Data  editing  functions 


Below  we  describe  the  algorithms  that  were  implemented 
to  edit  the  maps  prior  to  storing  them  as  quadtrees.  These 
algorithms  were  implemented  in  the  programming  language  C  to 
run  on  the  PDP-11/45,  VAX  11/780,  and  GRINNELL  configura¬ 
tion,  which  is  described  in  the  Appendix  to  this  report. 
The  maps  provided  were  initially  hand-digitized  and  stored 
in  a  chain  code  representation,  as  described  in  Section  3.1. 
Although  quadtrees  could  have  been  built  directly  from  the 
chain  code  representation,  it  was  considered  useful  instead 
to  use  picture  files  as  an  intermediate  representation 
between  the  initial  chain  codes  and  the  final  quadtrees. 
This  allowed  access  to  many  standard  routines  that  are  part 
of  our  software  library.  Below,  we  describe  the  nonstandard 
routines  that  were  used  on  this  project. 

The  algorithm  descriptions  proceed  in  the  following 
manner.  First  we  describe  the  program  MKBIN,  which  converts 
the  chain  codes  to  picture  file  format.  Then  three  routines 
for  manipulating  the  picture  files  are  described.  These  are 
FIXPIX  (changes  the  value  of  pixels  referenced  by  their 
coordinates) ,  RELABEL  (translates  one  list  of  pixel  values 
into  another) ,  and  LINERM  (erases  lines  from  maps) . 

The  function  MKBIN  (make  binary  array)  takes  as  input  a 
file  that  describes  a  chain  code  segmentation  of  the  origi¬ 
nal  maps  and  creates  a  picture  file.  For  the  following 
descriptions,  it  is  simplest  to  view  a  picture  file  as  a 
binary  array  that  f.as  been  laid  out  on  a  disk  in  a  row  by 
row  manner.  The  file  that  results  from  MKBIN  will  describe 
the  map  as  a  white  map  broken  up  by  a  series  of  black  lines. 
Since  the  black  lines  are  described  by  the  chain  codes, 
MKBIN  simply  traces  the  chain  code  on  a  binary  array,  mark¬ 
ing  each  pixel  that  lies  on  the  chain  code  as  black. 

Having  created  the  picture  file,  two  minor  utility  rou¬ 
tines  were  found  useful.  One  of  these  is  FIXPIX,  which 
changes  the  value  of  a  pixel  when  given  the  coordinates  of 
the  pixel  and  the  new  value.  This  is  the  equivalent  of 
assigning  to  an  entry  in  a  2-d  array.  FIXPIX  is  used  to  fix 
problems  that  result  from  errors  in  the  entry  of  the  chain 
code  digitization  and  also  errors  that  result  from  labeling 
the  regions  of  the  map.  The  other  utility  is  RELABEL,  which 
produces  a  copy  of  the  input  picture  file  where  all  the 
pixel  values  are  changed  according  to  a  given  translation 
table.  Both  of  these  utilities  are  used  for  changing  pixel 
values,  but  FIXPIX  does  this  based  on  specified  coordinates, 
whereas,  RELABEL  changes  all  pixels  that  have  a  given  value. 

The  land-use  classes  (as  well  as  contours,  etc.)  are 
labeled  by  using  a  standard  connected  component  program  and 
then  using  RELABEL  to  merge  the  labels  of  components  of  the 
same  class.  There  still  remains  the  problem  of  which  labels 


to  assign  the  pixels  that  lie  on  the  boundaries  of  the 
regions.  It  is  this  assignment  plus  the  original  hand  encod¬ 
ing  of  the  map  that  accounts  for  any  errors  in  the  calcula¬ 
tion  of  statistics  by  the  quadtree  algorithms.  The  assign¬ 
ment  decision  is  rather  arbitrary  (but  applied  consistently 
to  each  pixel  alike)  and  implemented  by  the  function  LINERM. 
The  decision  of  which  region  to  assign  a  given  boundary 
pixel  is  based  on  the  direction  of  easiest  movement  through 
a  picture  file  (which  is  from  right  to  left  and  from  top  to 
bottom) .  Thus  a  BLACK  pixel  (a  pixel  having  the  color  of 
the  boundary  line)  is  given  the  value  of  its  neighbor  on  the 
right  if  that  neighbor  is  not  BLACK  and  the  value  of  the 
neighbor  above  otherwise.  Note  that  at  any  point  during 
LINERM,  the  algorithm  stores  two  rows  of  the  picture  in 
core.  When  working  with  the  first  row  of  the  picture,  the 
second  neighbor  used  is  the  neighbor  below  since  there  is  no 
neighbor  from  above.  This  processing  is  repeated  over  the 
entire  picture  file  as  many  times  as  the  maximum  of  the 
line's  thickness  (measured  in  pixels) .  A  line  of  thickness 
greater  than  one  can  be  created  by  MKBIN  when  many  lines 
touch. 


Algorithm  j.1.  FIXP1X 
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/*  Take  a  picture  and  a  data  file  "input"  which  contains  records 
declared:  "x-coord  y-coord  newval" .  The  pixel  at  coord 
<x-coord  ,  y-coord >  will  be  changed  to  newval .  The  records  must 
be  in  ascending  order  of  y-coords  as  only  one  pass  is  made 
through  the  picture  (getting  new  rows  as  necessary)  . 

The  procedure  getrow  reads  tne  picture  file  filling  the  buffer 
with  the  next  row.  */ 

tixpix ( inpic ,  input,  numcols  ,  numrows) 

INTEGER  numcols,  numrows; 

DATA  FILE  inpic,  input; 

{ 

INTEGER  ARRAY  rowbuf f[ numcols} ; 

INTEGER  rownum  *  0; 

INTEGER  x  ,  y ,  val ; 

getrow  (inpic  ,rowbuf  f )  ; 

WHILE (NOT  end  of  file  "input") 

i 

getre cord  (input  ,x  ,y  ,val) ; 
if(y  >  rownum) 

FOR (r own um=r own um  TO  y) 
getrow (inpic , rowbuf f) ; 
rowbuf  f[x]  *  val; 

i 
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/ilgoritftin  3.2.  LInERR 


I*  Remove  the  black  pixels  from  a  multi-color  picture .  For  every 
pixel  in  the  picture,  do  the  following:  If  the  pixel  is  black, 
then  if  the  right  neighbor  is  not  black,  give  the  pixel  the 
value  of  its  right  neighbor.  If  the  right  neighbor  is  also 
black,  then  give  the  pixel  the  value  of  the  neighbor  above  it. 
Some  pixels  may  remain  black  (they  have  both  neighbors  black)  , 
if  so  the  algorithm  should  be  repeated.  */ 

linerm(inpic  .numcols  .numrows ) 

INTEGER  numcols  .numrows; 

DATA  FILE  inpic; 

{ 

INTEGER  ARRAY  inbuf fC 23C numcols+1 3 ; 

INTEGER  ARRAY  outbuffC numcols 3 ; 

INTEGER  POINTER  currpnt  ,  otherpnt; 

INTEGER  i  ,j  ; 

currpnt  =  0; 
otherpnt  =  1 ; 
getrow(inpic  , inbuf f[ 03 ) ; 
getrow(inpic  , inbuf fC 03) ; 
reset  inpic  file  to  beginning; 

/*  As  a  special  case  ,  the  top  row  actually  uses  the  neighbor  below 
it  rather  than  tne  neighbor  above  it  ,  and  the  right  hand  col 
uses  the  neighbor  to  the  left.  This  is  done  by  putting  an 
imaginary  row  above  the  first ,  and  an  extra  col  to  the  right 
of  the  last .  */ 

inbuf  ft  03t  numcols  3  “  inbuf  ft  03C  numcols-23  ; 
for(i=l  TO  numrows) 

( 

currpnt  =  (currpnt  *»  0) ; 

otherpnt  =  (otherpnt  ■*  0) ;  /*  flip  these  two  pointers  */ 

/*  Currpnt  always  points  to  the  current  row.  Otherpnt  points  at 
the  -row  above  .  */ 

getrow( inpic  ,inbuf  f[  currpnt 3 ) ; 

inbuf  f[currpnt3C numcols  3  ■  inbuf  f[  currpnt 3C  numcols-23  ; 

FOR(j*0  TO  numcols-1) 

t 

IF(inbuftCcurrpnt3[  j3  sa*  BLACK) 

IF(inbuf f[ currpnt3C j+13  <>  BLACK) 
outbuffC j3  ■  inbuf f[ currpnt3l j+13 ; 

ELSE 

outbuffC j3  *  inbuf fC other pnt 3l j3 ; 

ELSE 

outbuffC  j3  ■  inbuf  fC  currpnt3C  j3  * 
output (outbuff) ; 

} 


Algoritrun  3.3.  b*\£>lb* 
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/*  fciake  a  binary  array  from  a  set  of  chaincocUss .  The  chaincode 
used  is  as  follows:  " <coord-part><directional-part>»" 
<coord-part>  is  simply  a  six  digit  number  with  the  3  digit 
x-coord  and  the  3-digit  y-coord . 

<directional-part>  is  one  or  more  occurrences  of: 

<direction-character >  or  "C  <number><direction-character>]"  . 
<number>  is  a  2-digit  number  which  means  thac  the 
<direction-character >  occurs  number  times. 

<direction-character>  is  one  of  of: 

i  o  p 

\  1/ 

k  -*-  ; 

/l\ 

.  •  / 


If  the  character  is  ”k"  this  would  indicate  one  step  west, 
"  would  indicate  one  step  southwest,  etc.  These  symbols 
were  chosen  because  of  their  location  on  the  keyboard.  */ 


make bin (width  height) 

INTEGER  width,  height; 

/*  Create  a  binary  array  of  size  width  X  height.  */ 

/*  Tne  function  getchar  returns  the  next  character  that  is  not  a 
line-feed  from  file  "input" .  */ 

{ 

blwAKY  ARRAY  arr[ width ][ he ight ] ; 

INTEGER  xcoord  ,  ycoord  ,  numb  ,  i ; 

CHARACTER  ch; 


WHILE ( get coord s( xcoord , ycoord) ) 

{  /*  for  each  chaincode  in  the  file ...  */ 

arrC xcoord]C ycoord]  ■  1; 
cn  *  getchar (input) ; 

WHILE (ch  <>  '  t‘  ) 

( 

I*’(ch  ==  *C  '  ) 
run  (numb  ,ch) ; 

ELSE 

numb  *  1; 

FOR(i=l  TO  numb) 

CASE  OF  ch 

{ 


'  i '  : 

xcoord 

m 

xcoord 

- 

1 

ycoord 

X 

ycoord 

- 

l; 

'  o'  : 

ycoord 

s 

ycoord 

- 

1 

'  P'  * 

xcoord 

m 

xcoord 

+ 

1 

ycoord 

X 

ycoord 

- 

1; 

'  k'  : 

xcoord 

m 

xcoord 

- 

1 

*  •  *  • 
f  • 

xcoord 

m 

xcoord 

+ 

1 

•  •  • 

9  • 

xcoord 

m 

xcoord 

- 

1 

ycoord 

X 

ycoord 

+ 

1; 

•  • 

ycoord 

m 

ycoord 

+ 

1 

’  /'  * 

xcoord 

X 

xcoord 

♦ 

1 

ycoord 

* 

ycoord 

+ 

1; 

} 

arr[ xcoord][ ycoord]  *  1; 
cn  ■  getcnar (input) ; 

} 

) 

) 


BOOLEAN  FUNCTION  get coords (x ,y ) 

INTEGER  X ,  y/ 

/*  It  the  input  file  is  empty  return  FALSE.  Otherwise  read  the 
cooras  from  the  input  file  (into  x  ,y)  and  return  TRUE.  */ 

PROCEDURE  myget (numb  ,ch) 

INTEGER  numb; 

CHARACTER  ch; 

/*  Read  tne  2-digit  number  and  the  following  character  from  the 
input  file  ,  then  skip  the  character  ,  returning  the  number 
in  numb  and  the  character  in  ch.  */ 


Algorithm  3.4.  RELABEL 
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/*  Cnange  the  value  of  the  pixels  in  a  picture  as  determined  by  the 
labels  given  in  file  "labels"  .  This  file  has  as  its  first  value 
an  integer  which  is  the  largest  value  occurring  in  the  original 
picture,  followed  by  records  of  the  form  "old-val  new-val"  .  * / 

re  label  (inpic  ;iabels  .numcols  inumrows) 

INTEGER  numrows  ,  numcols; 

INTEGER  ARRAY  inpi ctnumcols]t numrows 3  ; 

DATA  FILE  labels; 

{ 

INTEGER  val  |new  ,old  ,i  ;  j  ; 

INTEGER  POINTER  table; 

val  =  getnum(labels) ; 

table  =  create-storage( ( val+1 )  *  sizeof  val); 

/*  Create -storage  is  a  system  function  which  dynamically  reserves 
the  number  of  words  given  by  the  parameter.  The  sizeof  operator 
returns  the  number  of  words  used  by  the  variable  following.  */ 

/*  Initialize  table  so  that  the  new  label  will  be  the  same  as  the 
old  label  ,  unless  a  change  is  indicated  in  the  file  "labels" .  */ 
FOR(i  *  0  TO  val) 

tablet i]  *  i; 

WHILE (not  at  end  of  "labels") 

{ 

old  *  getnum( labels ) ; 
new  s*  getnurn(!a  jcIs)  ; 
tablet  old]  *■  .*ew; 

} 

/*  Change  picture .  */ 

FOR(i=0  TO  numcols) 

FOR(j=0  TO  numrows) 

inpic[i]Cj]  *  tablet  inpi cti]tj]3: 
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4^.  Quadtree  encoding 
4.1.  Introduction 


This  section  describes  the  quadtree  encoding  algorithms 
as  well  as  various  primitive  functions  used  in  conjunction 
with  quadtree  data  structures.  Display  of  quadtree-encoded 
data  was  particularly  facilitated  by  the  ability  of  the 
GRINNELL  to  accept  specifications  of  rectangles  to  be  out¬ 
put.  It  should  be  noted  that  all  of  the  following  algor¬ 
ithms  work  on  the  digitized  version  of  the  maps  described  in 
Section  3  and  that  no  new  errors  are  introduced  by  these 
algorithms"  manipulations  of  the  quadtrees ,  since  the 
representation  of  the  digital  data  remains  exact.  No  devia¬ 
tion  from  pure  quadtree  representation  has  been  introduced. 

The  algorithm  descriptions  proceed  in  the  following 
manner.  First  we  present  a  set  of  primitive  functions  that 
form  the  building  blocks  for  later  algorithms.  Then  we  dis¬ 
cuss  two  algorithms  that  were  instrumental  in  building  the 
quadtree  database  from  the  digitized  maps.  The  first  of 
these  two  algorithms  builds  a  quadtree  from  a  map  by  scan¬ 
ning  the  map  in  a  row  by  row  fashion  (referred  to  as  raster 
scanning) .  The  second  algorithm  labels  the  connected  com¬ 
ponents  of  a  map. 

4.2.  Primitive  functions 

The  functions  SON,  FATHER,  SONTYPE,  NODETYPE,  BLACK, 
WHITE,  and  GRAY  can  be  thought  of  as  defining  the  quadtree 
as  an  abstract  data  type.  Although  their  implementation  is 
trivial,  their  usage  gives  the  other  quadtree  algorithms  a 
certain  independence  from  the  chosen  representation  of  the 
quadtree  data  structure.  Since  it  is  our  intent  to  experi¬ 
ment  with  other  quadtree  representations,  this  will  save 
future  programming  effort.  Currently  each  node  of  the  quad¬ 
tree  is  represented  by  a  record  consisting  of  five  pointers 
and  an  integer.  The  pointers  are  used  to  link  to  other 
nodes;  one  pointer  links  to  the  node's  father  and  the 
remaining  four  pointers  link  to  the  node's  four  sons  and  are 
indexed  by  the  quadrant  in  which  the  son  lies.  A  value  of 
NIL  is  stored  to  indicate  the  absence  of  a  son  in  a  given 
direction.  An  integer  value  is  used  to  uniquely  identify  the 
polygon,  land-use  class,  or  contour  to  which  the  region 
represented  by  the  node  belongs.  If  this  value  is  not  unique 
for  the  region,  then  the  value  is  considered  gray  (this  term 
comes  from  the  usage  of  gray  nodes  in  black  and  white 
binary-valued  quadtrees) . 

Using  such  a  quadtree  representation,  the  above  defin¬ 
ing  functions  work  as  follows.  The  function  SON  takes  a 
node  and  a  quadrant  as  parameters  and  returns  the  node  that 
is  the  son  of  the  given  node  in  the  given  quadrant  by  deref¬ 
erencing  the  appropriate  pointer.  Similarly,  the  function 


FATHER  takes  a  node  as  parameter  and  returns  the  father  of 
the  given  node  by  simply  dereferencing  the  appropriate 
pointer.  The  function  SONTYPE  takes  a  node  as  parameter  and 
returns  the  quadrant  that  expresses  the  direction  from  the 
father  of  the  given  node  to  the  given  node  by  comparing  the 
address  of  the  given  node  to  the  address  of  each  of  its 
father's  sons.  This  function  returns  a  special  value  NIL  to 
indicate  that  the  given  node  is  the  root  of  a  quadtree  and 
hence  has  no  father.  The  function  NODETYPE  takes  a  node  as 
its  parameter  and  returns  the  integer  data  item  that  is 
stored  at  that  node  which  generally  indicates  a  region 
color,  class  type,  or  elevation.  The  predicates  BLACK, 
WHITE,  and  GRAY  each  take  a  node  as  parameter  and  return 
true  if  the  value  of  NODETYPE  is  to  be  interpreted  as  having 
the  value  indicated  by  the  function's  name.  This  allows  mul¬ 
ticolor  quadtrees  to  be  easily  interpreted  as  binary-colored 
quadtrees  when  it  is  convenient  to  do  so. 

The  functions  OPSIDE,  CCSIDE,  ADJ ,  REFLECT,  QUAD,  and 
OPQUAD  provide  a  simple  set  of  operations  to  manipulate 
directions.  There  are  two  important  classes  of  directions 
used  by  quadtree  algorithms.  The  first  is  the  four  basic 
directions  denoted  N,  E,  S,  and  W  that  are  used  to  indicate 
the  side  of  the  square  that  lies  in  that  direction  from  the 
square's  center.  The  second  is  the  four  compound  directions 
denoted  NW,  NE,  SE,  and  SW  that  are  used  to  indicate  the 
quadrant  of  the  square  that  lies  in  that  direction  from  the 
square's  center.  The  functions  OPSIDE  and  CCSIDE  each  take 
a  side  as  parameter  and  return  respectively  the  side  in  the 
opposite  direction  and  the  side  in  the  direction  90  degrees 
counterclockwise  from  the  square's  center.  The  predicate 
ADJ  takes  a  side  and  a  quadrant  as  parameters  and  returns 
true  iff  the  given  quadrant  is  adjacent  to  the  given  side. 
For  example,  the  NE  quadrant  is  adjacent  to  both  the  N  and  E 
sides  but  not  to  the  S  or  W  sides.  The  function  REFLECT 
takes  a  side  and  a  qu_irant  as  parameters  and  returns  the 
quadrant  that  is  the  reflection  of  the  given  quadrant  with 
respect  to  a  line  through  the  center  of  the  square  that  is 
parallel  to  the  given  side.  For  example,  the  SW  quadrant  is 
the  reflection  of  the  NW  quadrant  with  respect  to  a  line 
through  the  square's  center  that  is  parallel  to  either  the  N 
or  S  sides.  The  function  QUAD  takes  as  parameters  two  sides 
and  returns  the  quadrant  that  is  adjacent  to  both  sides  if 
this  condition  uniquely  determines  one  quadrant.  If  it  does 
not  (i.e.,  the  two  sides  are  either  opposite  or  the  same), 
then  the  value  NEG  is  returned.  The  function  OPQUAD  takes  a 
quadrant  as  parameter  and  returns  the  quadrant  that  lies  in 
the  opposite  direction  from  the  center  of  the  square  (i.e., 
180  degrees) .  In  our  particular  implementation,  each  of  the 
two  classes  of  directions  is  represented  by  the  integers  0 
thru  3  inclusive;  so  the  above  functions  are  implemented  by 
modular  arithmetic  where  convenient  and  otherwise  by 
enumeration  of  the  possible  values  (i.e.,  table  lookup  via  a 
case  statement) . 


Central  to  the  approach  to  quadtrees  that  we  have 
adopted  is  the  ability  to  find  a  node's  neighbor  without 
storing  explicit  links  to  each  node's  neighbors  as  is  done 
in  some  other  implementations.  This  ability  is  encoded  in 
the  function  FIND_NEIGHBOR,  which  takes  as  parameters  a  node 
and  a  side  and  returns  a  node  that  abuts  the  indicated  side 
of  the  given  node  and  is  either  a  leaf  or  is  of  the  same 
depth  as  the  given  node.  The  manner  in  which  this  is  done  is 
described  in  the  tutorial  section  of  this  report.  The  func¬ 
tion  MAKEJNEIGHBOR  behaves  in  the  same  manner  as 
FIND_NEIGHBOR  except  that  if  it  fails  to  find  a  common 
ancestor  or  runs  into  a  leaf  before  it  has  finished  the  mir¬ 
rored  path,  then  it  modifies  the  tree  by  inserting  the 
sought-after  node  and  continues  on. 

The  remaining  functions  GETNODE,  CREATENODE,  and 
RETURNTOAVAIL  are  used  for  storage  management.  Unused  nodes 
are  kept  on  an  AVAIL  list.  The  function  GETNODE  returns  a 
used  node  by  first  looking  on  the  AVAIL  list  and  if  the 
AVAIL  list  is  empty  then  requesting  more  storage  from  the 
operating  system.  An  error  message  results  if  no  storage  is 
available  and  the  program  terminates.  The  function 
CREATENODE  takes  a  node,  a  quadrant,  and  an  integer  nodetype 
as  parameters  and  uses  GETNODE  to  create  a  new  node  of  the 
given  nodetype  which  has  the  given  node  as  father,  lies  in 
the  given  quadrant  of  the  given  node,  and  itself  has  no 
sons.  The  function  RETURNTOAVAIL  takes  a  node  as  parameter 
and  inserts  it  into  the  AVAIL  list. 

^.3.  Database  building 

Prior  to  constructing  the  quadtree  database,  the  maps 
were  stored  in  picture  files,  which  can  be  viewed  as  2-d 
arrays  laid  out  on  a  disk  in  row-by-row  order.  Thus  the 
first  task  to  be  performed  to  convert  each  picture  file  into 
a  quadtree  file,  which  is  a  preorder  listing  of  the  nodes  in 
a  quadtree.  This  is  accomplished  by  the  R2Q  (raster  to  quad¬ 
tree)  function.  This  function  reads  a  picture  file  one  row 
at  a  time  (raster  scan  order)  and  builds  the  corresponding 
quadtree  using  the  MAKE_NEIGHBOR  primitive.  As  the  quadtree 
is  being  built,  identical  leaf  brothers  are  merged  as  indi¬ 
cated  in  the  discussion  of  the  WINDOW  function.  An  addi¬ 
tional  efficiency  results  from  realizing  that  it  is  only 
necessary  to  check  for  these  mergers  on  even-  numbered  rows; 
any  leaf  on  an  odd-numbered  row,  still  has  two  brothers  that 
have  yet  to  be  read  in. 

The  original  picture  files  had  each  pixel  labeled 
according  to  the  land-use  class  (contour,  etc.)  to  which  it 
belonged.  Thus,  these  labels  were  the  only  distinctions 
that  could  be  carried  over  in  the  construction  of  the  quad¬ 
trees  by  R2Q.  However,  the  database  design  called  for  unique 
labels  on  each  connected  component  of  each  class.  Hence,  it 
was  necessary  to  perform  a  connected  component  analysis  in 


order  to  label  the  quadtrees  in  the  desired  format.  This  was 
done  by  the  function  QCONCOM  (quadtree  connected  component 
finder) .  This  analysis  was  performed  on  the  quadtree  data 
structure  directly  (instead  of  being  done  on  the  picture 
files  prior  to  quadtree  construction)  because  the  number  of 
nodes  to  be  processed  in  the  quadtrees  was  substantially 
smaller  than  the  number  of  pixels  to  be  processed  in  the 
original  picture  thereby  allowing  the  analysis  to  be  per¬ 
formed  faster.  The  function  QCONCOM  works  in  the  following 
manner  (processing  only  one  class  at  a  time).  The  first  step 
assigns  an  initial  tentative  labeling  to  the  quadtree.  This 
labeling  is  based  on  a  preorder  traversal  of  the  quadtree 
that  starts  in  the  northwest  corner  of  the  image  and  moves 
in  the  south  and  east  directions.  If  a  BLACK  node  is  met 
that  is  unlabeled  (with  respect  to  the  component  within 
which  it  is  contained) ,  then  a  new  label  is  created  for  it. 
When  processing  a  BLACK  node,  FIND  NEIGHBOR  is  used  to  exam¬ 
ine  its  southern  and  eastern  neighEors  to  determine  if  they 
are  also  BLACK,  but  have  no  component  label.  In  such  a  case, 
they  are  assigned  the  component  label  of  the  BLACK  node 
being  processed.  If  they  already  had  a  component  label,  then 
both  labels  are  placed  (as  an  ordered  pair)  on  an 
equivalence  list.  Once  all  the  nodes  have  been  tentatively 
labeled,  one  merges  the  equivalence  classes  and  then  updates 
the  component  labels  so  that  each  connected  component  has 
just  one  label. 


Algorithm  4.1.  PRIMITIVES 
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/*  The  following  is  a  description  of  the  primitive  functions 
used  in  the  quadtree  algorithms.  */ 

node  F UWCT ION  s  on ( p  ,  i ) 

/*  Given  node  p  and  quadrant  i  ,  return  the  node  which  is  the 
son  representing  quadrant  i  of  node  p.  */ 

node  FUNCTION  father (p) 

/*  Given  node  p,  return  the  node  which  is  the  father  of  p.  */ 

INTEGER  FUNCTION  sontype(p) 

/*  Given  node  p,  return  q  where  son(father (p)  ;q)  «  p.  If  p  is 
the  root,  then  return  NIL.  */ 

INTEGER  FUNCTION  nodetype (p) 

/*  Return  the  value  of  node  p.  This  can  be  considered  as  GRAY, 
WHITE;  or  BLACK  for  a  binary  tree;  and  GRAY,  WHITE  or  a  class 
type  or  elevation  level  value  for  a  multi-colored  tree  .  */ 

bOOLEAN  FUNCTION  black (p) 

/*  TRUE  when  nodetype (p)  is  BLACK  if  the  tree  is  binary;  or 

when  nodetype(p)  is  a  value  specified  as  BLACK  by  the  user 

it  the  tree  is  multi-colored.  */ 

BOOLEAN  FUNCTION  White (p) 

/*  TRUE  when  nodetype(p)  is  WHITE  if  the  tree  is  binary,  or 

wnen  nodetype(p)  is  a  value  specified  as  WHITE  by  the  user 

if  the  tree  is  multi-colored.  */ 

BOOLEAN  FUNCTION  gray(p) 

/*  TRUE  iff  nodetype (p)  is  GRAY.  */ 

INTEGER  FUNCTION  opside(b) 

/*  Given  a  side  b,  return  the  opposite  side  (e.g.; 
opside(E)  ■  W) .  */ 

INTEGER  FUNCTION  ccside(b) 

/*  Returns  the  side  adjacent  to  side  b  in  the  clockwise 
direction  (e.g.,  ccside(E)  *  N) .  */ 

BOOLEAN  FUNCTION  adj(b,i) 

/*  TnUE  iff  quadrant  i  is  adjacent  to  boundary  b  of  the  node's 
block  (e.g.,  adj(N,Nw)  ®  TRUE;  adj(N,Sw)  **  FALSE).  */ 

INTEGER  FUNCTION  reflect(b,i) 

/*  Returns  the  quadrant  which  is  adjacent  to  quadrant  i  along 
boundary  b  (e.g.,  ref  le  ct  (N  ,Nw)  -  Sw)  .  */ 

INTEGER  FUNCTION  quad(b,c) 

/*  Returns  tne  quadrant  bounded  by  b  and  c  it  it  exists  and  the 
value  NEG  if  it  does  not  exist.  */ 

INTEGER  FUNCTION  opquad(q) 

/*  Returns  the  quadrant  opposite  (non-ad jacent )  to  q 
(e.g.,  opquad(NW)  ■  SW)  .  */ 
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noae  FUNCTION  find_neigmx>r  (q  ,s) 

J*  Return  the  node  which  is  adjacent  to  side  "s"  or  node  "q" 
and  is  either  a  lear  or  is  at  the  same  depth  as  node  "q" . 

This  is  done  by  following  the  father  links  until  the  common 
ancestor  is  reached  and  then  following  the  rerlected  path 
downward,  stopping  s,»ort  only  if  a  leaf  is  met.  */ 

node  POINTER  q; 

INTEGER  s; 

t 

node  POINTER  p; 

INTEGER  i  .stypeq; 

/*  first  fina  a  common  ancestor.  */ 

IE (NULL (son type (q) ) )  /*  Common  ancestor  does  not  exist.  */ 

RETURN  ( NIL }  ; 

ELSE  IF(ad j (s  ,sontype(q) ) )  /*  Neighbor  is  not  a  sibling  - 

go  up  to  next  level .  * / 
p  «■  find__neighbor  ( father  (q)  ,s); 

ELSE  /*  Neighbor  is  a  sibling,  father  is  a  common  ancestor.  */ 
p  »  father (q); 

/*  After  finding  the  common  ancestor  ,  reflect  about  side  "s" 
back  to  the  level  of  tne  original  request.  */ 

IF (NULL (p)  OR  NULL( son( p  .reflect ( s  ,  son type (q) ))) ) 

/*  Either  there  was  no  common  ancestor  or  p  is  a  leaf 
and  in  either  case  p  is  what  we  want  to  know; 
so  ,  don' t  change  it .  */ 

RETURN (p) ; 

ELSE  /*  Return  the  calculated  son.  */ 

RETURN ( son (p  )reflect(s  ,sontype(q) ) ) ) ; 

i 

node  FUNCTION  make_neighbor (q  ,s ) 

/*  Return  the  node  which  is  adjacent  to  side  "s"  of  node  "q" 
and  at  the  same  depth  as  the  node  "q"  .  Tnis  is  done  by 
following  the  path  through  the  tree  that  would  lead  us 
to  said  neighbor  it  it  existed  and  creating,  along  the  way, 
any  nodes  tnat  are  necessary,  whenever  such  nodes  are  created  , 
all  created  sons  are  set  to  WHITE.  They  are  later  reset  to 
GRAY  or  BLACK  as  appropriate,  c.f.,  find^ neighbor  */ 

node  POINTER  q; 

INTEGER  s; 

l 

node  POINTER  p; 

INTEGER  i  .stypeq; 

/*  First  find  the  nearest  common  ancestor  .  */ 

IF(NULL(sontype(q) ) )  /*  Common  ancestor  does  not  exist.  */ 

t 

/*  Create  a  common  ancestor  and  initialize  its 
pointers .  */ 

p  •  createnode (NULL  , NULL  .GRAY)  ? 
stypeq  ■  quad ( ccside ( s )  .opside ( s ) ) ? 
p->sonsC stypeq]  ■  q; 


q->fathr  »  py 

/*  Create  the  other  three  eons  of  p.  */ 
createnode (p  ,opquad(  stypeq)  ,wHITE)y 
createnode ( p  ,opquad( re fleet  Is  , stypeq) )  .whlTE) y 
createnode (p  ,reflect(s  .stypeq)  iWHITE) y 

ELSE  IF (ad j (s  , son type (q) ) )  /*  Neighbor  is  not  a  sibling  - 

go  up  to  the  next  level.  */ 
p  *  make_neighbor (father (q) ,s)y 
ELSE  /*  Neighbor  is  a  sibling.  Father  is  common  ancestor.  * 
p  *  father (q); 

/*  After  finding  the  nearest  common  ancestor  ,  reflect  about 
side  "s"  back  to  the  level  of  the  original  request.  */ 

IF ( NULL ( son (p  ,ref  lect(s  ,sontype(q) ) ) ) ) 

i  /*  If  the  node  does  not  have  children  to  descend  a  level  , 
change  the  node  to  gray  and  give  it  children.  */ 
p->  node  type  =  GRAY ; 

FOR  ( i  «  NW  ,NE  ,SE  ,SW ) 
createnode  ( p  ,i  .nHITE)  ; 
return(son(p  ,reflect(s  ,sontype(q) ) ) ) ; 


node  FUNCTION  getnode() 

/*  Reserves  storage  for  a  quadtree  node  and  returns  a  pointer 
to  this  unit  of  storage  */ 

node  FUNCTION  createnode (root  ,s  ,t) 

/*  Create  a  node  p  with  nodetype  t  which  corresponds  to  son  s 
of  node  root  and  return  p.  * / 
node  POINTER  root ; 

INTEGER  s  ,t ; 

l 

node  POINTER  p; 
p  *»  getnode  ( )  y 
if (root  I*  NIL) 

root->8onsCs}  *  py 
p->fathr  -  rooty 
p->ntype  *  ty 
for  (i  •  nw;ne,se;sw) 
p->sonsCi]  *  NILy 
return(p)  y 

} 

PROCEDURE  returntoavail(p) 

/*  Return  node  p  to  the  available  storage  pool.  */ 


Algontrua  <*.2.  QCOwiCOE 
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/*  Run  a  connected  components  algorithm  on  a  binary  quadtree  - 
i  .e  assign  every  connected  component  a  unique  label.  This  is 
achieved  in  three  steps .  Step  1  assigns  labels  to  each  BLACK 
node  .  This  is  done  by  traversing  the  tree  and  tor  every  BLACK 
node,  examining  its  eastern  and  southern  neighbors.  If  these  are 
unlabeled,  a  new  label  is  generated  for  the  current  node.  If 
either  of  tnem  are  labeled  ,  and  the  current  node  is  unlabeled  , 
then  assign  current  node  the  label  of  its  neighbor.  It  the 
neighbor  is  labeled  and  the  current  node  is  labeled  ,  then  these 
labels  are  equivalent  and  so  the  pair  of  labels  is  added  to  a 
list  of  equivalence  classes .  The  second  step  is  to  put  the  list 
of  equivelance  classes  into  a  hierarchal  order  so  that  all  of 
the  nodes  in  the  class  can  be  given  one  label.  No  algorithm  is 
given  for  this  step  -  for  an  example  of  a  typical  algorithm  of 
tnis  Kind  see  Knuth ,  Vol  1 .  The  third  step  simply  traverses 
the  tree  again,  relabeling  each  node  to  the  value  of  the 
representative  for  its  equivalence  class .  * / 

component (quadtree ) 

/*  "Quadtree"  is  a  pointer  to  the  input  tree.  At  the  end  of  this 
algorithm,  "quadtree"  will  point  to  the  labeled  tree.  */ 
node  POINTER  quadtree ; 

1 

pairlist  POINTER  merges;  /*  Pointer  to  the  list  of  pairs  of 

equivelances  .  * / 

merges  =  NIL; 

labe 1( quadtree ) ;  /*  step  1  */ 

Process  the  equivalences  in  the  list  merges;  /*  step  2  */ 
update ( quadtree ) ;  /*  step  3  */ 


label  ( p) 

/*  Perform  step  1.  Assigns  labels  to  node  p  and  its  sons.  */ 
node  POINTER  p; 

{ 

node  POINTER  q; 

INTEGER  i; 

If (gray(p) ) 

FOR(  i  *  nw,ne;se,sw) 
labe  1( son (p  ,i) ) ; 

ELSE  If ( black (p) 

{ 

q  •  find_ neighbor (p  ,E) ; 
if ( NOT ( NULL ( q ) ) ) 

label_ad  jacent(q  ,NW  ;sw  ;p) ; 
q  *  find_neighbor  (p  ;s) ; 
if ( NOT ( NULL ( q ) ) ) 

labelad  jacent  (q  ,NW  ,NE  ,P) ; 
if (NOT (labe led (p) ) ) 

p->nodetype  *  getnewregion( ) ; 
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lately  ad  jacent(r  ;ql  ,q2  ,p) 

/*  Find  all  descendants  of  node  r  adjacent  to  node  p  -  i .e .  , 
in  quadrants  ql  and  q2 .  */ 
node  POINTER  r  ,p; 

INTEGER  ql  ,q2; 

i 

IF(gray(r ) ) 

{ 

label_ad  jacent(son(r  ,ql )  ,ql  ,q2  ,p) ; 
label' ad  jacent( son (r  ,q2)  ,ql  ,q2  ,p) ; 

} 

ELSE  IF(black(r ) ) 

assign  label ( p  ,r ) ; 


as3iqn_label  (p  ,q) 

/*  Assign  a  label  to  nodes  p  and  q  if  they  do  not  already  have  one  . 

If  both  have  different  labels  ;  then  enter  them  in  “merges"  .  */ 
node  POINTER  p  ',q ; 

l 

IF(labeled(p)  AND  labeled (q)) 

l 

IF (node type (p)  <>  nodetype (q)) 

add  <nodetype(p)  ,nodetype(q)  >  to  merges; 

} 

ELSE  IF(labeled(p) ) 

q-> node  type  *  nodetype(p); 

ELSE  IF ( labeled (q)) 

p->node cype  »  nodetype (q); 

ELSE 

p->nox_a  type  ■  q->  node  type  »  getnewregion( )  y 

i 


update (p) ; 

/*  Perform  step  3  .  */ 
node  POINTER  p; 

l 

INTEGER  i; 

IF(gray( p) ) 

FOR(i  «  nw;ne;se;sw) 
update  ( son  ( p  ;i ) ) ; 

ELSE  IF(black(p) ) 

p->nodetype  »  equivalence  of  value  nodetype (p); 


) 


Algor  it  run  4.3.  R2U 
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/*  Convert  an  input  picture  in  the  form  of  a  binary  array  (or 

raster)  into  a  binary  quadtree.  Basically;  the  algorithm  works 
by  doing  a  raster  scan  of  the  input  picture  ,  and  as  each  pixel 
is  read  ;  the  quadtree  is  modified  so  that  it  would  be  a  valid 
quadtree  representing  the  input  picture  if  all  unprocessed  pixel 
were  WHITE.  This  is  in  contrast  to  an  algorithm  which  first 
builds  a  complete  quadtree  with  one  node  per  pixel  and  then 
attempts  to  merge  nodes  (replace  GRAY  nodes  that  have  all  sons 
the  same  color  with  a  node  of  that  color).  Tne  input  picture  is 
read  one  row  at  a  time  by  a  special  function  getrow  which 
returns  the  next  row  of  the  picture .  The  function  color  returns 
the  color  of  the  pixel  given  as  an  argument.  The  boolean 
function  lastrow  is  true  iff  the  current  row  is  the  last  row  of 
tne  picture,  nmenever  all  the  children  of  a  node  have  been 
processed  ;  an  attempt  is  maae  to  merge  them  togetner  .  because 
or  this  there  is  a  distinction  between  odd  rows  and  even  rows 
(no  pixels  in  an  odd  row  can  ever  complete  tne  processing  of  all 
tne  children  of  a  gray  node  ;  hence  tnere  is  never  an  attempt  to 
merge  after  processing  any  of  tnese  pixels).  The  picture  is 
assumed  to  be  an  2**N  by  2**N  picture  -  if  not;  WHITE  pixels  are 
assumed  to  fill  it  out.  */ 

node  POINTER  quadtree ( p  .width ) ; 

/*  Grven  a  picture  p  (viewed  as  a  list  of  rows)  and  its  width, 
return  a  quadtree.  */ 

LIST  p; 

INTEGER  width; 

{ 

BOOLEAN  ARRAY  qCl:width];  /*  Holds  a  row  of  the  picture.  */ 
node  POINTER  first; 

INTEGER  i; 

q  »  getrow(p); 

first  *  createnode (NIL  ,  NULL  ,q[l] ) ;  /*  First  pixel.  */ 
oddrow(q  .first  .width) ; 
i  «  2; 

p  *  NEXT ( p) ; 

first  *  evenrow(getrow(p)  .make  neighbor  ( first  ,S)  .i.widtn); 

WHILE  (NOT  lastrowO)  DO 

{  /*  Process  the  rest  of  the  rows.  */ 

p  -  NEXT ( p) ; 

oddrow (getrow (p)  ifirst  .width); 
p  -  NEXT(p); 
i  »  i+2; 

^  first  *  evenrow(ge  trow(p)  ,make_neighbor  ( first  ,S)  ,i  .width); 
while (NOT  NULL ( father ( first) ) ) 

first  »  father ( first) ;  /*  Set  first  to  root  of  the  tree.  */ 

re turn ( first ) ; 

1 


oddrow  (row  ,nd  .width) 

/*  Add  the  odd-numbered  row  of  width  "width"  represented  by  array 
"row"  to  a  quadtree  whose  node  "nd"  corresponds  to  the  first 
pixel  in  the  row.  */ 
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INTEGER  width; 

INTEGER  ARRAY  rowC 1 : width] ; 
node  POINTER  nd; 

l 

nd->nodetype  =  color (rowCl] ) y 
F0R(i=2  UNTIL  width) 

l 

nd  ■  make_neighbor  (nd  ,E) ; 
no->nodetype  *  color (row[ i ]) ; 

) 


node  FUNCTION  evenrow(row  , first  ,i  fwidth) 

/*  Add  even  numbered  row  "i"  of  width  "widtn"  represented  by  array 
"row'*  to  a  quadtree  whose  node  "first"  corresponds  to  the  first 
pixel  in  tne  row.  During  this  process,  merges  of  nodes  having 
four  sons  of  the  same  color  are  performed.  */ 

INTEGER  ARRAY  row; 
node  POINTER  first; 

INTEGER  i  .width; 

( 

node  POINTER  p)r; 

INTEGER  j; 

p  «*  first; 

IF ( NOT  lastrowO) 

/*  Remember  the  first  node  of  the  next  row.  */ 
first  *  make  neighbor (p  ;S) ; 

FOR(j-l  UNTIL  width-1 ) 

r  »  makej^neighbor  (p  ;e)  ; 
p-> node type  »  color (rowC j] ) ; 

IF (EVEN ( j ) ) 

merge ( i  ; j  .father ( p) ) ; 

P  -  r; 

) 

p->nodetype  =  color  (rowC  width] )  ;  /*  Don't  invoke  make_^neighbor  lor 

the  last  pixel  in  a  row.  */ 

IF (EVEN (width) ) 

merge (i  )width  ;father(p) ) ; 

RETURN ( first ) ;  /*  Return  the  first  node  of  the  next  row.  */ 

) 


node  FUNCTION  merge(i;j)p) 

/*  Attempt  to  merge  a  node  having  four  sons  of  the  same  color 
starting  with  node  "p"  at  row  "i"  column  "j".  */ 
node  POINTER  p; 

INTEGER  i  ,j ; 

{ 

INTEGER  k; 

WHILE ( EVEN ( i )  AND  EVEN ( j )  and 

(nodetype ( 8on(p  ,NW) )  ■  node type( son ( p  ,NE) )  ■ 
node t ype (son (p  ,SE) )  ■  node t ype ( son (p  ;sw )) ) 
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Tabulation  of  results 


Below  we  describe  the  tables  of  data  collected  about 
the  quadtree-represented  regions.  The  times  are  measured  in 
seconds  by  a  special  routine  on  the  VAX  11/780,  which  allows 
us  to  factor  out  the  disk  I/O  time.  Thus  the  times  reported 
reflect  the  notion  of  CPU  time  *s  implemented  on  that 
machine.  These  times  were  measured  while  the  machine  was 
not  loaded  with  any  other  jobs,  because  the  timing  routine 
does  vary  in  its  results  as  the  system  load  changes.  For 
many  algorithms,  a  better  idea  of  the  cost  in  time  can  be 
determined  from  such  machine  independent  concepts  as  number 
of  nodes  visited.  These  are  also  included  in  many  of  the 
tables  or  are  easily  deducible  from  the  descriptions  of  the 
algorithms  in  the  algorithm  overview. 

The  tables  are  discussed  in  the  order  that  they  appear 
appended  to  this  section.  These  tables  treat  the  data  base 
as  a  collection  of  possibly  unconnected  regions  on  the  maps 
that  are  logically  connected  by  the  sharing  of  some  pro¬ 
perty,  e.g.,  having  the  same  land-use  class. 

The  first  group  of  tables  (4.1-3)  are  the  QUADTREE 
BUILDING  STATISTICS.  These  reflect  a  process  by  which  the 
R2Q  algorithm  was  used  to  build  a  separate  quadtree  for  each 
logically  connected  class  of  polygons  is  a  picture  file.  No 
execution  time  is  indicated  because  the  times  were  dominated 
by  the  cost  of  reading  an  entire  picture  file,  one  line  at  a 
time.  Each  quadtree  took  approximately  3  minutes  to  build. 
If  the  quadtree  had  been  constructed  by  building  a  complete 
4-ary  tree  and  then  merging  where  possible,  it  would  be 
necessary  for  the  memory  to  be  large  enough  to  contain 
262,144  (512x512)  nodes.  By  merging  during  the  building 
process, as  R2Q  does,  a  much  smaller  maximum  memory  is 
required  in  practice.  The  size  required  is  recorded  in  the 
column  'nodes  created' .  Note  that  never  are  more  than 
15,000  nodes  needed.  The  following  column  indicates  the  per¬ 
centage  of  this  maximum  that  was  actually  used  when  the  tree 
was  finished.  The  remaining  columns  give  a  breakdown  of 
the  number  of  nodes  of  each  type  in  the  resulting  binary 
quadtree. 

The  CONNECTED  COMPONENT  RESULTS  (Tables  4.4-6)  record 
the  data  collected  on  three  variations  of  the  connected  com¬ 
ponent  algorithm,  QCONCOM.  in  each  variation,  the  algorithm 
uses  two  of  the  neighbors  of  a  node,  as  described  in  the 
algorithms  overview,  to  assign  a  tentative  label  to  a  node. 
The  'number  of  neighbors  sought'  is  the  number  of  times  this 
process  of  finding  a  neighbor  must  be  performed,  thus  yield¬ 
ing  an  indication  of  its  importance  to  the  algorithm's 
analysis.  The  three  variations  are  three  different  methods 
of  finding  the  required  neighbors.  For  purposes  of  com¬ 
parison,  the  average  cost  for  a  single  finding  of  a  neighbor 
is  calculated  to  show  clearly  the  variance  within  each 
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technique  as  well  as  the  relative  costs  among  techniques. 
The  average  cost  is  measured  in  number  of  nodes  accessed  per 
neighbor  found,  since  the  ammount  of  work  performed  by  the 
algorithm  is  proportional  to  the  number  of  nodes  accessed, 
this  average  cost  measure  gives  an  accurate  view  of  the 
relative  tradeoffs  among  the  various  methods.  The  portion 
of  a  second  required  to  find  a  neighbor  (an  alternative 
measure)  could  not  be  calculated  due  to  the  inaccuracy  of 
the  system  timing  algorithm.  Also,  measurement  in  seconds 
of  algorithm  efficency  can  be  misleading  because  the  algor¬ 
ithms  were  coded  in  a  highlevel  language  and  some  of  the 
timing  differences  could  reflect  the  relative  efficencies  of 
the  compiler's  optimizer  rather  than  that  of  the  quadtree 
algorithm. 

The  first  method  is  FINDNBR,  which  is  the  FIND_NEIGHBOR 
primitive  mentioned  in  the  algorithm  overview  and  described 
in  the  tutorial  section.  The  time  in  the  final  column  is 
for  this  method.  The  second  method,  ROPES,  is  due  to 
Hunter's  quadtree  work  referred  to  in  the  tutorial.  It  con¬ 
sists  of  placing  a  link  directly  between  each  neighbor  of 
the  same  size.  This  results  in  a  reduction  in  execution 
time  at  a  major  cost  in  storage  (due  to  storing  the  extra 
links)  .  The  third  method  was  discovered  during  work  on  this 
project  and  consists  of  causing  the  traversal  algorithm  to 
pass  as  parameters  the  neighbors  of  each  subtree's  root. 
This  requires  more  time  than  ropes,  but  does  not  require  as 
much  memory.  The  added  memory  cost  with  respect  to  the 
F IND__NE  IGHBOR  technique  results  from  the  additional  stack 
size  needed  due  to  the  larger  parameter  list  of  the  recur¬ 
sive  routines.  The  average  value  is  based  on  equating  the 
cost  of  passing  a  parameter  to  a  subroutine  with  the  cost  of 
dereferencing  a  link.  This  equivalence  is  clearly 
compiler-dependent  as  well  as  machine-dependent.  The  final 
column  of  the  table  indicates  overall  execution  time  of  the 
connected  component  algorithm,  QCONCOM,  which  uses  the 
FIND  NEIGHBOR. 
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5.  Region  analysis  and  manipulation 
5.1.  Region  analysis 

The  functions  described  in  this  section  are  used  to 
gather  basic  statistics  about  the  regions  encoded  by  the 
quadtree  data  structure.  The  simplest  of  these  is  NDCOUNT, 
which  sets  three  global  variables  to  indicate  the  number  of 
gray  nodes,  white  nodes,  and  black  nodes  in  the  given  quad¬ 
tree.  This  is  achieved  by  performing  a  preorder  traversal 
of  the  quadtree  —  i.e.,  first  it  calculates  its  statistic 
for  the  current  node  (in  this  case  incrementing  a 
global  counter)  and  then  it  recursively  processes  the 
node's  four  subtrees  (if  they  exist) .  It  should  be  noted 
that  the  functions  described  in  this  and  in  subsequent 
sections  are  currently  implemented  as  stand  alone  pro¬ 
grams  that  are  invoked  by  executing  a  file  under  the  UNIX 
operating  system.  This  entails  a  fairly  large  amount  of 
housekeeping  details,  such  as  decoding  the  arguments  with 
which  the  file  is  invoked,  opening  various  files,  and 
initializing  various  devices.  None  of  this  will  be  discussed 
herein,  nor  will  it  appear  in  the  algorithm  descrip¬ 
tions.  Among  the  other  details  that  are  thus  swept  under 
the  rug,  so  to  speak,  would  be  the  initializing  of 
the  variables  queried  by  the  functions  BLACK,  WHITE,  and 
GRAY  in  order  to  determine  (when  processing  a  multicolored 
quadtree)  which  nodes  should  be  considered  of  the  indi¬ 
cated  colors.  Although  we  speak  herein  of  functions 
computing  values,  we  actually  have  programs  that  generate 
files  and  output  listings  containing  function  values. 
For  instance,  the  implementation  of  NDCOUNT  terminates 
by  outputting  the  counts  for  the  three  types  of  nodes. 

Closely  related  to  NDCOUNT  is  a  function  called  AREA. 
AREA  takes  as  parameters  a  node  and  an  indicated  width  for 
that  node.  AREA  calculates  the  area  and  centroid  of  the 
region  encoded  by  black  nodes  relative  to  this  indicated 
size  for  the  entire  quadtree.  This  is  achieved  by  a  preorder 
traversal  that  works  as  follows.  If  the  current  node  is  a 
leaf,  then  its  size  is  added  to  the  global  count  in  accor¬ 
dance  to  whether  or  not  it  is  black.  If  the  current  node  is 
not  a  leaf,  then  AREA  processes  each  of  its  four  subtrees 
using  a  width  value  adjusted  to  half  of  the  value  associated 
with  the  current  node. 

Next  in  order  of  complexity  is  HANDW ,  which  takes  as 
its  parameters:  a  node,  the  x  and  y  coordinates  of  its  upper 
left  corner,  and  its  width,  its  value  is  the  coordinates  of 
the  upper  left  corner,  the  height,  and  the  width  of  the 
smallest  rectilinear  rectangle  (i.e.,  the  smallest  rectangle 
with  sides  parallel  to  the  x  and  y  axis)  that  encloses  all 
the  regions  that  are  considered  black.  This  is  done  by  com¬ 
paring  the  coordinates  of  each  black  node  to  the  most 
extreme  values  found  so  far.  Again  we  are  dealing  with  a 


preorder  traversal  of  the  quadtree  with  the  x,  y,  and  width 
values  being  updated  as  one  descends  from  a  gray  node  to  its 
children. 

The  last  of  the  statistics  gathering  functions  is  PER¬ 
IMETER.  It  also  updates  a  global  variable  and  calculates 
its  desired  value  via  a  preorder  traversal.  Like  AREA,  it 
takes  a  node  and  its  width  as  parameters.  Unlike  AREA,  it 
uses  the  width  value  to  calculate  the  length  of  a  node's 
side  instead  of  the  node's  area.  The  function  PERIMETER 
returns  as  its  value  the  sum  of  the  lengths  of  the  perime¬ 
ters  of  all  the  black  regions.  Like  AREA,  it  does  this  from 
the  point  of  view,  so  to  speak,  of  the  black  nodes.  The 
side  of  a  black  node  is  part  of  the  perimeter  (and  hence  its 
length  is  to  be  counted)  only  if  the  neighbor  on  that  side 
of  the  black  node  is  a  white  node.  The  neighbor  is  located 
using  FIND  NEIGHBOR. 


algorithm  b  .1 .  KbCGUistf 


96 


/*  A  simple  tree  traversal  to  count  the  number  of  GRAi'  ,  wril'ft.  ana 
BLACK  nodes  .  * / 

INTEGER  numgray  -  0; 

INTEGER  numwhite  «  0; 

INTEGER  numb lack  »  0; 

PROCEDURE  ndcount(rt) 
node  POINTER  rt; 

{ 

INTEGER  i; 

IF(gray(rt->ntype ) ) 

l 

numgray  »  numgray  +1; 

FOR(  i  ■  WW  ,NE  ,SE  ,SW) 

ndcount (rt->sons[ i] ) ; 

i 

ELSE 

IF(black(rt->ntype) ) 

numb  lack  **  numblack  +  1 ; 

ELSE 

numwhite  »  numwhite  +  1; 
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/*  Given  a  quadtree  ,  compute  the  area  and  centroid  of  the  black 
region  of  that  tree .  */ 

INTEGER  area,  xsum  ,  ysum; 


main (root ,  width) 

/*  To  compute  the  area  ,  simply  sum  up  the  number  of  pixels  in  each 
black  node  . 

To  compute  the  centroid  ,  for  x-coord  sum  up  all  of  the  x-coord 
of  black  pixels  and  divide  by  area;  for  y-coord  sum  y-coords 
and  divide  by  area .  */ 
node  POINTER  root ; 

INTEGER  width; 

l 

area  *  xsum  *  ysum  *  0; 

doarea ( root , width  ,0  ,0) ;  /*  compute  area 

(stored  in  global  variable  area)  */ 
xcent  »  xsum/area;  /*  xcoord  of  centroid  */ 

ycent  *  ysum/ area;  /*  ycoord  of  centroid  */ 


doarea (root ,  width,  f x  ,  fy) 

/*  This  function  does  the  work  of  computing  the  area  and  the 

centroid.  For  each  black  node,  add  The  number  of  pixels  to  the 
global  variable  area;  sum  up  the  x-coordinates  and  add  to  the 
global  variable  xsum ,  and  sum  up  the  y-coordinates  and  add  this 
to  global  variable  ysum.  */ 
node  POINTER  root; 

INTEGER  width  ,  temp; 

INTEGER  f x ,  fy;  /*  Coords  of  the  upper  left  pixel  of  the  node.  */ 
if (gray(root) ) 

l  /*  for  each  child  ,  compute  area  */ 

doarea  (width/ 2  ,son(root  ,NW)  ,fx  ,fy) ; 
doarea  (width/ 2  ,son(root  ,NE)  ,fx+width/2  ,fy) ; 
doarea (width/ 2  ,son(root  ,SE)  ,fx+width/2  ,fy+width/2) ; 

^doarea (width/ 2  ,son(root  ,SW)  ,fx ,fy+width/2) ; 

else 

if (black (root) ) 

1 

/*  Incorporate  area  of  current  node  into 
running  totals .  */ 
temp  »  width  *  width; 

xsum  *  xsum  +  (fx  +  width/2  -  .5)  *  temp; 
ysun  »  ysun  ♦  (fy  width/2  -  .5)  *  temp; 
area  ■  area  +  temp; 
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'  Find  the  smallest  enclosing  rectangle  for  the  black,  area  of  a 
quadtree.  Specify  this  rectangle  by  its  upper  left  coordinates 
and  its  height  and  width.  */ 

ITEGER  leastx  .leasty  .highx  .highy;  /*  The  highest  and  lowest 

values  yet  found.  */ 


tin  (root  .width) 

1  Given  a  quadtree  ,  call  handw  to  get  the  highest  and  lowest 
values  of  x  and  y  coords.  The  upper  left  corner  is  <least-x  , 
least <*y>  and  the  width  and  height  is  the  difference  between  the 
x' s  and  y' s  respectively.  */ 
lode  POINTER  root; 

LftTEGEK  height  .width; 

Leastx  *  leasty  »  width  +  1; 
lighx  *  highy  =*  0; 
iandw(root  ,0  ,0  .width) ; 
leight  *  highy  -  leasty; 

/idth  *  highx  -  leastx; 


indw(root  ,x  ,y  .width) 

*  X  and  y  are  the  coordinates  of  the  upper  left  corner  of  the 
node  width  is  the  width  of  the  node  .  If  the  node  is  black , 
check  if  the  extreme  corners  of  the  node  are  within  the  least 
and  high  bounds  -  if  not,  then  change  the  bounds.  If  the  node 
is  gray,  check  the  sons.  */ 

INTEGER  x, y, width ; 
lode  POINTER  root; 


Lf (gray (root ) ) 

{  /*  For  each  son,  do  handw.  */ 

handw  (son  (root  ,Nw)  ,x  ,y  .width/ 2) ; 
handw  (son  (root  ,NE)  ,x+width/2  ,y  .width/ 2) ; 
handw (son (root  ,SE)  ,x+width/2  ,y+width/2  .width/ 2) ; 
handw (son (root  ,SW)  ,x  ,y+width/2  .width/ 2) ; 

) 

ilse 

if (black (root) ) 

l 

if(x  <  leastx)  leastx  ■*  x; 

if ( (x+width-1)  >  highx)  highx  -  x  +  width  -  1; 
if(y  <  leasty)  leasty  -  y; 

if ( (y+width-1)  >  highy)  highy  -  y  +  width  -  1; 

) 
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/*  Given  a  quadtree  and  its  width ,  compute  the  length  of  the 

perimeter  of  the  black  areas .  This  is  done  by  traversing  the 
tree  and  calling  addperim  for  each  black  node .  Addperim  looks 
at  each  of  the  neighbors  of  the  black  node.  If  that  neighbor  is 
white  ,  then  the  length  of  the  edge  is  added  to  the  perimeter 
total .  If  it  is  gray ,  then  sons  along  the  edge  of  the  the  black 
node  are  run  with  addperim.  */ 

INTEGER  perimlength  *  0; 
perim(root  .width) 

/*  Traverse  the  tree  calling  addperim  for  black  nodes.  */ 
node  POINTER  root; 

INTEGER  width; 

l 

If (gray (root ) ) 

FOR(i  =  a W  to  sw) 

perim(width/2  .son (root  ,i) ) ; 

ELSE 

If (black (root) ) 

l 

addperim ( find  neighbor  (root  ,N)  .width  ,SE  ,Sw) ; 
addperim ( f ind_neighbor  (root  ,E)  .width  ,NW  ,SW) ; 
addperim  ( find_  neighbor  (root  ,S)  .width  ,NW  ,NE) ; 
addperim ( f ind_neighbor  (root  ,w)  .width  ,NE  ,SE)  ; 

} 

addperim  (  root  .width  ,ql  ,q2) 

/*  Root  is  a  neighbor  of  a  black  node  .  If  root  is  white  .  add  width 
to  the  perimeter  length,  if  it  is  gray;  then  perform  addperim  on 
the  children  which  are  adjacent  to  the  original  black  node 
(quadrants  ql  and  q2) .  */ 
node  POINTER  root; 

INTEGER  width  ,ql  ,q2 ; 


If (nil(root) )  /*  The  black  node  was  on  the  edge  of  the  tree  -  no 

neighbor  exists .  */ 
perimlength  =  perimlength  +  width; 

ELSE 

IF (white (root) ) 

perimlength  ■  perimlength  4  width; 

ELSE  If (gray(root) ) 

{ 

addperim  (son  (root  ;ql)  .width  ,ql  ,q2) ; 
addperim  (son  (root  ,q2)  .width  ;ql  ,q2) ; 

1 
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5.2.  Region  manipulation 


Of  the  basic  operations  described  herein,  the  only  one 
that  does  not  produce  a  new  quadtree  is  the  PT2P0LY  func¬ 
tion,  which  given  a  quadtree,  a  coordinate  structure  (x  and 
y  coordinates  of  the  upper  left  corner  and  the  width  for  the 
entire  tree) ,  and  a  (u,v)  coordinate  pair,  returns  the  value 
(color)  of  the  leaf  node  that  represents  the  region  that 
contains  the  coordinate  pair.  Unlike  the  statistics  gather¬ 
ing  programs  that  had  to  traverse  the  entire  tree,  the 
PT2P0LY  function  only  visits  those  nodes  that  lie  on  a 
direct  path  between  the  quadtree's  root  and  the  sought  after 
leaf.  This  is  done  by  determining  the  quadrant  of  the 
current  node  within  which  the  coordinate  pair  lies  and  then 
recursing  down  into  that  subtree  while  updating  the  coordi¬ 
nate  structure  to  reflect  the  new  location. 

One  of  the  basic  operations  that  take  one  quadtree  as  a 
parameter  and  return  a  new  quadtree  as  the  result  is  the 
WINDOW  function.  Besides  its  quadtree  parameter,  the  WINDOW 
function  also  takes  a  specification  of  where  the  window 
should  be  placed—  i.e.,  the  current  width  of  the  quadtree, 
the  coordinates  of  the  upper  left  corner  of  the  window,  and 
the  width  of  the  window.  For  the  present,  the  window  must  be 
a  square  whose  width  is  a  power  of  2.  The  new  quadtree  is 
constructed  by  recursively  performing  the  following  steps. 
Find  the  smallest  subtree  of  the  given  quadtree  that  con¬ 
tains  the  window.  If  this  subtree  coincides  with  the  window 
then  return  the  subtree.  If  this  subtree  is  a  leaf  then 
return  a  leaf  of  the  same  color.  Otherwise,  split  the 
current  window  into  quadrants  and  process  each  of  these 
subwindows  with  respect  to  the  current  subtree,  upon  return¬ 
ing  from  each  recursive  call,  it  is  necessary  to  check  if 
four  leafs  are  brothers  of  the  same  color,  and  when  this 
happens,  replace  the  father  by  one  of  the  four  leafs.  This 
process  results  in  a  quadtree  that  represents  the  windowed 
portion  of  the  map  encoded  by  the  given  quadtree. 

The  two  basic  operations  that  take  two  quadtrees  as 
parameters  are  the  set-theoretic  operations  of  INTERSECTION 
and  UNION.  Both  of  these  functions  work  on  binary  quadtrees, 
taking  two  quadtrees  as  parameters  and  creating  a  resulting 
quadtree.  In  both  cases,  we  assume  that  the  input  quadtrees 
are  of  the  same  width  and  have  the  same  upper  left  coordi¬ 
nates.  If  this  were  not  the  case,  the  user  could  perform 
the  WINDOW  function  to  align  the  two  quadtrees.  Like  the 
statistics  gathering  functions,  these  two  operations  perform 
preorder  traversals  of  the  quadtree  parameters.  However, 
now  the  traversals  are  performed  in  parallel?  so  that  at  any 
time  during  the  processing,  the  algorithms  keep  track  of  the 
two  nodes  (one  in  each  quadtree)  which  correspond  to  the 
same  areas  in  the  two  encoded  maps.  The  logic  of  these  two 
operations  is  summarized  below. 


if  current 
nodes  are 


function  is 
INTERSECTION 


function  is 
UNION 


both  black 
both  white 
both  gray 


one  black 


one  white 


return  black 

return  white 

recurse 

return  other 
subtree 

return  white 


return  black 
return  white 


recurse 


return  black 


return  other 
subtree 


In  the  above  table,  "recur se"  indicates  that  one  needs  to 
traverse  each  of  the  remaining  subtrees,  and  "return  other 
subtree"  indicates  that  the  value  of  the  function  is  a  copy 
of  the  other  subtree.  Just  as  with  windowing,  when  the 
recursion  unwinds,  one  has  to  check  to  see  if  four  brothers 
are  identical  leaves  and  merge  them  as  indicated  in  the 
discussion  of  the  WINDOW  function.  Examples  of  performing 
these  basic  set-theoretic  operations  are  shown  in  Figures 
5. 1-5. 3.  Table  5.1  shows  the  area  of  the  landuse  polygons 
in  Figure  5.3,  using  the  naming  conventions  discussed  in 
Section  5.3. 

The  set-theoretic  operation  of  complement  can  be  per¬ 
formed  using  the  more  general  function  QMASK.  QMASK  takes  a 
quadtree  and  a  range  as  parameters  and  returns  a  quadtree 
that  has  all  the  nodes  with  values  within  the  range  set  to 
BLACK  and  all  the  nodes  with  values  outside  the  range  set  to 
WHITE.  Thus  the  tree  that  results  from  QMASK  is  always  a 
binary-colored  quadtree.  The  QMASK  algorithm  is  implemented 
as  a  preorder  traversal  of  the  input  quadtree  that  simul¬ 
taneously  constructs  the  output  quadtree.  In  constructing 
the  output  quadtree,  the  QMASK  algorithm  merges  nodes  when 
necessary  as  indicated  in  the  discussion  of  the  WINDOW  func¬ 
tion. 

The  final  quadtree  manipulation  function  is  QDISPLAY, 
which  does  double  duty  both  as  a  quadtree  manipulator  (it 
truncates  quadtrees)  and  as  an  output  routine.  The  parame¬ 
ters  of  QDISPLAY  are  a  quadtree,  specifications  of  how  the 
quadtree  should  be  displayed  (location,  width,  coloring 
algorithm,  etc.),  and  the  depth  at  which  the  quadtree 
should  be  truncated.  The  coloring  algorithm  is  determined 
by  two  flags,  COLOR  and  BLOCK.  If  COLOR  is  true,  then  the 
quadtree  is  displayed  as  is,  each  node"s  value  being  inter¬ 
preted  as  a  color.  If  COLOR  is  true  and  BLOCK  is  false, 
then  the  quadtree  is  output  as  a  binary-colored  quad¬ 
tree,  where  the  colors  mapped  to  BLACK  are  defined  by  set¬ 
ting  the  range  used  by  the  primitive  function  BLACK.  If 
BLOCK  is  true,  then  a  special  table  of  colors  is  used 


and  the  colot  of  a  node  is  determined  by  its  depth  and  its 
binary-colored  value.  If  COLOR  is  false,  then  one  has 
the  option  of  setting  the  maximum  depth  of  a  node  that  will 
be  displayed.  When  a  gray  (internal  non-leaf)  node  is 
to  be  displayed,  the  function  examines  the  gray  node's  des¬ 
cendants  and  considers  the  node  to  be  BLACK  if  the  black 
nodes  (when  weighted  according  to  their  depth  in  the  tree) 
exceed  the  white  nodes  (when  similarly  weighted) .  Thus  a 
node  is  output  as  BLACK  in  the  truncated  tree,  if  it  is 
BLACK  or  it  is  a  gray  node  at  the  maximum  depth  and  the 
average  color  of  the  region  it  subtends  is  more  black  than 
white.  The  algorithm  is  implemented  as  a  preorder 
traversal  of  the  input  quadtree  that  outputs  the  nodes  in 
the  order  they  are  visited.  Figures  5. 4-5. 8  show  the 
output  of  QDISPLAY  when  COLOR  and  BLOCKS  are  false  and  the 
polygon  flood. center  is  considered  BLACK.  These  figures 
give  an  idea  of  the  initial  gentle  degradation  of  the  image 
as  the  quadtree  is  truncated.  Table  5.12,  discussed  in 
Section  5.3,  shows  that  Figure  5.5  uses  only  two-thirds 
the  number  of  nodes  as  Figure  5.4,  with  virtually  no  loss 
in  the  basic  image  shape.  This  shows  that  quadtree  trunca¬ 
tion  is  a  useful  image  approximation  technique. 


.2 .  Kesult  of  executing  lwYi.a.^tx.TIun  on  the  entire 
land-use  map  and  tne  complement  of  tne  stn 
elevation  level  (<*U0-i>d0  ft.  elevation)  of  tne 
topography  map . 


if  Figure  5.3.  Result  of  executing  INTERSECTION  on  the  1st 

elevation  level  (0-100  ft.  elevation)  of  the 
topography  map,  the  flood. center  region  of  the 
floodplain  map,  and  the  entire  land-use  map. 
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Figure  5.6.  hesult  of  executing  gDltoPLAY  on  flood. center  of 
the  flood-plain  map  using  to  levels. 
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/*  Given  a  tree  ana  tne  coordinates  of  a  point ,  return  tne  value  of 
tne  leaf  node  containing  that  point.  It  is  assumed  that  in  a 
database  system  the  tree  used  for  this  operation  would  have  a 
different  node  value  for  each  polygon.  This  would  allow  a 
simple  lookup  to  determine  the  corresponding  polygon  once  the 
node  value  in  th-?  tree  has  been  found  by  this  algorithm.  */ 

INTEGER  FUNCTION  pt2poly(fx.  f y ,  width,  x coord  ,  ycoord  ,  rt) 

node  POINTER  rt;  /*  l*ne  root  of  the  tree.  */ 

INTEGER  f x  ,  f y;  /*  This  is  the  upper  left  coord  of  the  tree  .  */ 

INTEGER  width;  /*  The  width  of  the  tree.  */ 

INTEGER  xcoord ,  ycoord;  /*  The  given  coord  being  searched  for  */ 

IF(1  gray(nodetype (rt ) ) ) 

RETURN ( node type (rt ) ) ; 

ELSE  /*  Gray  tree  -  find  which  quadrant  contains  the  coord.  */ 

IF (ycoord  <  (fy  +  width/2))  /*  North  half  */ 

IF (xcoord  <  (ix  +  width/2))  /*  NW  */ 

RETURN (pt2poly(fx  ,fy  .width/ 2  .xcoord  .ycoord  ,son(rt  ,Nw) )  ; 

ELSE  /*  NE  */ 

RETURN (pt2poiy( f x+width/ 2  ,fy  ,width/2  .xcoord  iycoord  , 

8on(rt  ,NE ) )  ; 

ELSE  /*  South  half  */ 

IF(xcoord  >*  (fx  -r  width/2))  /*  SE  */ 

RETURN(pt2poly(fx+width/2  ;fy+width/2  ,widtn/2  .xcoord .ycoord  ; 

son(rt  ,SE) ) ; 

ELSE  /*  SW  */ 

RETURN (pt2poly(fx  ,fy+width/2  .width/ 2  .xcoord  .ycoord , 

1  son(rt.SW)); 


Algontnra  s.b.  window 
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/*  Given  a  quadtree  and  its  width  ;  along  with  the  speciiications  of 
a  window i  create  a  tree  which  is  the  section  of  tne  input  tree 
specified  by  the  window.  */ 

main  (root  iwidtn  ,wcol  ,wrow  iwwidth) 
node  POINTER  root;  /*  Root  is  the  input  tree.  */ 

INTEGER  width  ,wcol  iwrow  iwwidth; 

/*  width  is  the  width  of  the  input  tree  ,  <wcol  ;wrow>  is  the  coord 
of  the  upper  left  corner  of  the  window;  wwidth  is  the  size  or 
the  window  (must  be  a  power  of  2) .  */ 

l 

node  mode;  /*  Dummy  node  to  start  the  answer  tree.  */ 
node  POINTER  rptr;  /*  Eventual  root  of  the  answer  tree.  */ 

rptr  =  mode ; 
mode  .fathr  *  NIL; 

dowindow(root  iwidth  ,0  iO  .rptr  ,NW  .wwidth  ,wcol  iwrow) ; 
rptr  »  son  ( rptr  ,Nw) ; 

/*  Rptr  is  now  the  root  of  the  answer  tree  .  */ 

) 


dowindow( xnrt  .inwidth  iincol  , inrow  ioutfthr  iwhichson  ioutwidth  i 
out col  .outrow) 

/*  From  the  input  tree  create  a  tree  described  by  the  given  ’window 
which  has  father  outfthr  and  is  son  whichson .  */ 
node  POINTER  inrt  ioutf thr ;  /*  Inrt  is  the  root  of  the  current 
input  tree .  Outfthr  is  the  father  of  the  output  tree .  */ 
INTEGER  inwidth  iincol  .inrow  iwhichson  ioutwidtn  ioutcol  ioutrow; 

/*  Inwidth  i  incoi  and  inrow  are  tne  window  described  by  the  input 
tree ;  outwidth  i  outcol  i  and  outrow  are  the  window  to  be 
described  by  the  tree  being  built  and  whichson  indicates  the 
son type  of  tne  tree  being  built  in  relation  to  tne  whole 
output  tree .  */ 

i 

node  POINTER  nd; 

INTEGER  i  i  goodquad; 

/*  First  i  cut  the  input  tree  down  to  the  smallest  subtree  which 
contains  the  window  desired  -  i.e.i  if  the  window  is  completely 
contained  in  one  of  the  children  of  the  input  tree  i  make  the 
the  input  tree  that  child  (and  continue  for  its  children...)  */ 
goodquad  »  1 ; 

while  ( (goodquad  <>  NEG)  AND  (inrt  GRAY)) 

i 

goodquad  *  NEG; 

/*  For  each  quadrant  i  check  if  window  in  quadrant .  */ 

if  (inrect( incoi  iinrow  iinwidth/2  ioutcol  ioutrow  ioutwidth) ) 
goodquad  ■  NW; 

if  (inrect(incol+inwidth/2  iinrow  iinwidth/2  ioutcol  ioutrow  i 

outwidth ) ) 

goodquad  •  NE; 

if (inrect( incol+ inwidth/ 2  ,inrow+inwidth/2  .inwidth/ 2 .outcol ; 

outrow  ioutwidth) ) 

goodquad  «  SE; 

it  (inrect(  incoi  ,inrow+inwidth/2  .inwidth/ 2  .outcol  ioutrow  i 


outwidth ) ) 


good quad  *  SW; 
i  t  (  qoodq uad  <  >  NEG) 

l  /*  Window  in  input  quadrant  - 

make  input  tree  tnat  quadrant  */ 
inrt  *  son  ( inrt  ;goodquad)  ; 
inwidth  *  inwidth/ 2; 
if ( good quad  ==  NE  OK  S£) 
incoi  =  incol  +  width; 
if(goodquad  **  SE  OK  SW) 
inrow  »  inrow  +  width; 

} 

} 

if ((incol  “  outcol)  AND  (inrow  ==  outrow)  AwD 
(inwidth  outwidtn)) 

/*  If  the  windows  are  the  same  ;  then  the  output  tree  is  the 
same  as  the  input  tree  ,  so  copy  the  input  tree  .  */ 
copysub(inrt  .outfthr  jwhichson) ; 
else 

if (node type ( inrt)  <>  GRAY) 

/*  If  the  input  tree  is  a  leaf  node  ;  then  the  window  is 
also  a  leaf  of  the  same  color.  */ 
createnode  (outfthr  .whichson  inodetype  ( inrt ) ) ; 

else 

{  /*  No  single  child  of  the  input  tree  contains  the  window ;  and 
the  input  tree  is  not  a  leaf  node.  Therefore;  repeat 
dowindow  on  each  quadrant  of  the  windo"  desired.  To  do 
this  ;  first  install  a  GRAY  node  in  the  output  tree  and 
then  call  dowindow  for  each  quadrant.  */ 
nd  «  createnode  (outfthr  ;whichson  ;inrt->ntype) ; 
dowindow  (inrt  ;inwidth  ;incol  ,  inrow  ; 

nd  ;NW  ;outwidth/2  .outcol  loutrow)  ; 
dowindow  (inrt  .inwidth  iincol  ,  inrow  ; 

nd  ,NE  ;outwidth/2  ,outcol+outwidth/2  .outrow) ; 
dowindow  (inrt  ;inwidth  ;incol  ;  inrow  ; 

nd  ,SE  ;outwidth/2  ,outcol+outwidth/2  ,outrow+outwidth/2) 
dowindow  (inrt  .inwidth  ;incol  .inrow  ; 

nd  ;outwidth/2  .outcol  ;outrow+outwidth/2 ) ; 
if (All  children  of  nd  have  the  same  nodetype) 

t 

nd->ntype  ■  nodetype(son(nd  ;nw)  ) ; 

Return  all  children  of  nd  to  avail  list; 

) 

} 

} 

BOOLEAN  FUNCTION  inrect (bigcol  Ibigrow  .bigwidth  .litcol  .litrow  ; 

litwidth) ; 

/*  Return  TKUE  if  and  only  if  the  second  window  is  contained 
in  the  first  window.  */ 

INTEGER  bigcol  ;bigrow  .bigwidth  ;iitcol  .litrow  ;iitwidtn; 

if ((bigcol  <»  litcol)  AND  (bigrow  <■  litrow)  and 
(bigcol+bigwidth  >■  litcol+litwidth )  AND 
(bigrow-*- bigwidth  >»  litrow+litwidth) ) 
return(TKUL) ; 


} 


re turn (FALSE) ? 


copysub(root  ;fthr  ,whi chson ) 

/*  Create  a  copy  of  the  subtree  with  root  "root"  .  The  created 
tree  will  be  son  "whichson"  of  the  node  "fthr"  (part  of  the 
global  answer  tree )  .  */ 
node  POINTER  root, fthr? 

INTEGER  whichson; 

l 

node  POINTER  root ; 

INTEGER  i; 

nd  **  createnode  ( fthr  ;whi chson  .node type  (root )) ; 
for  ( i  -  NW  ,NE  ;se  ,Sw ) 

copysub(  son  (root  ii )  ,nd  ;i )  ? 


Algorithm  5.7.  INTERSECTION 
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/*  Given  two  (possibly  multi-colored)  quadtrees  ,  return  a  binary 
quadtree  which  is  the  intersection  of  the  input  trees .  */ 

node  FUNCTION  inter  (rtl  ,  rt2  ,  fthr  ,  whichson) 

/*  Return  the  intersection  of  the  trees  whose  roots  are  pointed  to 
rtl  and  rt2 .  This  is  done  by  simultaneously  traversing  tne 
input  trees  and  performing  inter  on  each  quadrant.  If  either 
tree  is  WHITE,  the  intersection  is  WHITE.  If  either  tree  is 
BLACK ,  the  intersection  is  the  other  tree . 

Fthr  is  the  father  of  the  resulting  tree  .  This  allows  the 
current  subtree  to  be  inserted  into  tne  complete  output  tree . 
whichson  tells  which  son  of  the  output  tree  the  current  tree 
is.  The  function  is  initially  called  with  these  values  equal 
to  NIL.  */ 

node  POINTER  rtl  ,  rt2 ,  fthr ; 

int  whichson; 

( 

node  POINTER  nd ; 

INTEGER  i; 

IF  (white (node type (rtl ) )  OR  white (nodetype (rt2) ) ) 

RETURN (createnode ( fthr  .whichson  .WHITE) ; 

IF  (black(nodetype(rtl ) ) 

RETURN (copysub(rt2  .fthr  .whichson) ) ; 

IF  ( black( nodetype (rt2) ) 

RETURN (copy s ub( rtl  .fthr  .whichson) ) ; 

/*  Both  trees  are  GRAY  -  call  inter  for  each  quadrant .  */ 

nd  »  createnode  (fthr  .whichson  .GRAY) ; 

FOK(i  -  NW.NE.SE.SW) 

inter  (rtl ->sons[i]  ,rt2->sons[  i]  ,nd  ,i); 

IF(all  children  of  nd  are  WHITE) 

{  /*  Condense  tree .  */ 

nd->ntype  ■  WHITE; 

FOR(i  -  NW  .NE  .BE  ,SW) 

returntoavail(nd->sons[i] ) ; 

} 

RETURN (nd); 


algorithm  b  .d  •  Uwlou 
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/*  Given  two  (possibly  multi-colored)  quadtrees  ,  return  a  binary 
quadtree  which  is  the  union  of  the  input  trees.  */ 

node  FUNCTION  union ( rtl  ;  rt2  ,  fthr  ,  whichson) 

/*  Return  the  union  of  the  trees  whose  roots  are  pointea  at  by 
rtl  and  rt2 .  Tnis  is  done  by  simultaneously  traversing  tne 
input  trees  and  performing  union  on  each  quadrant,  if  either 
tree  is  black,  tne  union  is  BLACK.  If  either  tree  is  WHITE, 
the  union  is  the  other  tree .  F thr  is  the  fatner  o,\  the 
resulting  tree .  This  allows  tne  current  subtree  tc  be  inserted 
into  the  complete  output  tree.  Whichson  indicates  tne  sontype 
of  the  current  tree  relative  to  the  output  tree.  Union  is 
initially  called  with  these  values  equal  to  NIL.  */ 

node  POINTER  rtl  ;  rt2  ;  fthr; 

int  whichson; 

{ 

node  POINTER  nd ; 

INTEGER  i; 

IP  (black(nodetype (rtl ) )  OR  black(nodetype(rt2) ) ) 

RETURN  (createnode(  fthr  iwhichson  ;BLACK } ; 

IF  (white(nodetype(rtl) )  /*  If  treel  is  WHITE,  */ 

RETURN(copysub(rt2  ifthr  ;whichson) ) ;  /*  copy  tree 2  .  */ 

IF  ( whi te( node type ( rt2 ) )  /*  If  tree2  is  WHITE,  */ 

RETURN  ( copysub(  rtl  ifthr  ^whichson) ) ;  /*  copy  treel.  */ 

/*  Both  trees  are  GRAY  and  union  must  be  applied  to 
each  quadrant .  */ 

nd  *  createnode (fthr  ;whichson  ;GRAY) ; 

FOR(i  »  ww;ne,se;sw) 

union ( rtl ->sons[ij  ;rt2->sonsCi]  ,nd  ;i) ; 

IF (All  children  of  no  are  BLACK) 
i  /*  Condense  tree  .  */ 
nd->ntype  ■  BLACK; 

FOR(i  »  nw,ne;se;sw) 

returntoavail(nd->sons[ i] ) ; 

) 

RETURN ( nd ) ; 


nlgoritnm  b.9.  givneiK 
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/ *  tor  a  quadtree  with  root  'root*  ,  for  eacn  leaf  node  ,  it  its 
value  is  between  tne  parameters  nigh  and  low  (inclusive)  ,  tnen 
set  node  value  to  BLACK  .  else  set  to  wril'i't .  unen  necessary, 
marge  cniloren  of  a  gray  none  toaetner  .  */ 

li .<  TfcObK  low  ,  hign ; 

qmasK(rt) 
nous  Puli.4i'tK  rt; 
l 

lftTbtibK  i  ; 

Ir(oKAi (rt) ) 
t 

to«(i  =  tjw  ,hb  ,t»u  ) 

qmasK(rt->sons[ i] )  ; 

IF (all  children  of  rt  are  tne  same  leaf  color) 
t  /*  condense  tree  */ 

rt->ntype  *  noaetype (rt->sons[nw] ) ;  •* 

FOK(i  =«  tnu  ,Bb  ,br.  ,bW) 

returntoavail (rt->sonsC i3 ) ; 
i 

} 

CiLbc, 

I*( (rt->ntype  >=  low)  AWli  (rt->ntype  <=*  high)) 
rt->ntype  »  BLACK; 
hLbh 

rt->ntype  *  wHIfb; 


main  (root  .low  ,high) 

node  FoIw'i'LK  root;  /*  input  tree  */ 

Iwi'bcbK  low,  high;  /*  input  parameters  -  values  in  this  range 

are  BLACK .  */  • 

{ 

qmasK(root ) ; 

Kb'i’UKK (root)  ; 


Algor  itnxu  5.1b.  ublbPLAY 
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/*  Lisplay  a  quadtree  on  the  grinnell  .  Tnis  is  done  by  traversing 
tne  tree  .  while  tne  traversal  is  being  ctone  .  tne  program  Keeps 
tracK  or  tne  size  and  position  or  tne  current  node  ,  and  displays 
that  node  on  tne  grinnell .  ‘i'nere  are  two  boolean  options  - 
color  and  olocK.  If  color  is  TkUu  ,  tnen  nodes  will  oe  displayed 
by  tneir  color  ,  ir  FALbL  tnen  all  non-wHlTb  nodes  are  displayed 
as  LLAbk.  Ir  block  is  true,  then  the  displayed  color  depends 
on  tne  depth  or  tne  node  in  tne  tree  -  not  its  value  .  t  most 
one  or  tnese  options  may  be  true.  Additionally,  it  color  is 
FaLpl  ,  tne  user  may  wish  to  display  nodes  only  down  to  a  certain 
deptn  by  adjusting  maxdepth.  In  tnis  case,  tne  smallest  node 
size  displayed  can  be  changed,  F or  example,  with  a  512  X  512 

picture  ,  there  are  10  levels  and  tne  smallest  node  is  one  pixel 

wide .  ir  the  user  sets  maxdepth  as  y  ,  then  tne  smallest  node  is 
2  a  2  pixels,  For  any  gray  nodes  at  level  y,  tne  runction  ok 
auus  up  tne  number  ot  black  pixels  ot  its  cnildren  and  ir  more 

tnan  nalr  are  black ,  tne  gray  node  is  displayed  as  black  ,  otner 

wise  it  is  displayed  as  white  . 

'next'  is  a  runction  which  returns  the  value  of  the  next  node 
of  tne  preoroer  traversal  of  the  input  tree  stored  in  tree .  */ 


DA1A  FILL  tree;  /*  preorder  traversal  ot  the  input  tree.  */ 
liaTtULA  color  ,  block; 
lbTLGLK  maxlevel; 

IdTLOLK  larrClO];  /*  Tnis  array  holds  grinnell  color  values  to  be 

used  with  option  block.  */ 


coloror(val  .curr level) 

/*  Letermine  tne  actual  color  value  to  be  displayed  on  tne  grinnell. 
for  tne  node  witn  value  val .  This  is  determined  by  the  options  , 
tne  value  of  the  node  and  possibly  tne  level  of  the  none  in  tne 
tree  .  */ 

IwTLGLk.  val  ,  currievel  /*  currlevel  is  the  level  in  the  tree  ot 

the  current  node .  */ 
t 

IhTLGEK  v; 

IF (block) 

{  /*  color  determined  by  level  */ 

v  «  larrC currlevel  -  1]; 
lF(black( val) ) 
v  »  v*25b; 

/*  wnite  nodes  will  be  displayed  as  some  tint  of  red  ,  black  nodes 
will  be  displayed  as  some  tint  of  blue  .  Multiplying  by  2:>o 
snifts  a  red  value  to  a  blue  one  when  displayed.  */ 
ret urn ( v) ; 

} 

bLbb 

IF (black (val ) ) 

IF (color) 

KbTUiUM  ( val ) ; 

LLbL 

KhTURw ( pLACa ) ; 
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KFTUkw  (  Miii'i'h )  ; 

i 


uK(wiatn) 

/*  iteturn  tne  number  or  pixels  of  tne  current  node  ana  its 
children  wnich  are  Dlack.  */ 
iuTnGtK  widt n; 
i 

XwTcidfiK  val; 

val  =  next();  /*  causing  this  function  to  have  a  side  efrect  */ 
It (blacK( val ) ) 

aL'i'uRlJ (width  *  width); 

It‘  (white  (val ) ) 

KtfUKW(O) ; 

/*  gray  node  -  calculate  black  pixels  in  chilaren  */ 

KiiTUtUji  (bk(width/2)  +  bk(widtn/2)  +  bk(width/2)  +  bk(  widtn/2 ) )  ; 

i 


display! rcol  ,rrow  ,width  .currlevel) 

/*  Display  tne  next  node  of  the  preorder  traversal  on  the  grinnell . 
This  node  has  its  upper  left  corner  at  fcol  ,  trow  ana  has  width 
widtn .  It  is  at  level  currlevel  in  the  tree.  */ 

.  IWTEvifcK  rcol,  frow,  width,  currlevel; 

l 

ItfTtbLK  col,  total,  val; 

* 

val  *  next! ) ; 

It'! gray! val)) 
t 

col  *  colorof! val  .currlevel ) ; 

<*«rite  to  grinnell  a  square  at  fcol  ,rrow  of  size  width  witn 
color  col.>  /*  this  is  a  command  and  not  a  comment  */ 

) 

fcLbfi 

IF(currlevel  ==  maxlevel) 

IF (block) 

{ 

col  *  colorof (val  .currlevel ) ; 

<write  to  grinnell  a  square  at  fcol  ,frow  of  size  width  and 
color  col . > 

1 

bL£>£ 

l 

total  »  Dk(width/2)+bk!width/2)+bk(wiath/2)+bk!width/2) ; 

1F( (2  *  total))  >  (width  *  width)) 
val  ■  1; 
fiLSF 

val  »  0; 

col  »  colorof (val  .currlevel ) ; 

<write  to  grinnell  a  square  at  fcol  ,frow  of  size  wiath  and 
color  col . > 


} 
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DLSt. 

I  /*  Display  children  or  a  gray  node  */ 
display! rcol  ,(frow+wiath/2)  ,(wioth/2)  ,  (currle vel+i )) ; 
display! ( lcol+widtn/2)  ,(frow+widtn/2)  ,(widtn/2)  ,(currlevel+i ) ) ; 
display! (fcol+width/2)  ,irow  ,(wioth/2)  ,lcurrlevel+l ) ) ; 
display! ! f  col  ,trow  , !  width/2)  ,! currlevel+i ) ) ; 

) 

i 


main!  f  col  ,rrow  , color  ,blocx  ,maxlevel  ,widtn)  ; 

iKl'cdhK  rcol  , trow  .width; 
i 

./*  Angle  brackets  enclose  commands  written  in  English  and  are  not 
just  comments.  */ 

<r’ill  larr  witn  grinnell-depenaent  values  used  witn  option  blocx .  > 
-display! rcol  ,trow  ,widtn  ,1 ) ; 
i 
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5.3.  Tabulation  of  results 

In  the  tables  presented  in  this  section,  the  basic 
unit  of  manipulation  is  the  connected  component.  The  names 
of  these  basic  units  are  created  by  suffixing  a  digit  to  the 
land-use  class  (or  contour)  to  which  the  unit  belongs.  These 
digits  can  be  dereferenced  by  referring  to  the  figures  in 
Section  3.  For  example,  in  Table  5.2,  the  first  polygon 
name  we  encounter  is  acc.l.  Looking  at  Figure  3.4,  one  sees 
the  19  components  of  the  class  ACC.  The  component  labeled  1 
in  that  figure  is  the  polygon  refered  to  by  the  name  acc.l. 
The  units  of  the  flood-plain  map  are  so  few  that  they  are 
given  names  of  their  own,  i.e.,  left  and  right  bank  instead 
of  bank.l  and  bank. 2.  Note  that  there  are  no  tables  showing 
the  execution  times  for  the  PT2P0LY  function.  This  is 
because  all  times  were  less  than  a  tenth  of  a  second  and 
hence  were  beyond  the  range  of  the  system  timing  algorithm. 

The  first  group  of  tables  (Tables  5.2-4)  are  the  AREA 
RESULTS  tables.  They  are  organized  according  to  which  map 
the  polygons  (i.e.,  simply-connected  components)  belong. 
They  summarize  the  results  of  two  programs:  NDCOUNT  (which 
counts  the  number  of  black  nodes,  i.e.,  those  belonging  to 
the  polygon)  and  AREA  (which  calculates  the  area  in  pixels 
and  centroid  (first  moment)  of  the  polygon) .  The  execution 
times  immediately  follow  the  results  of  the  same  algorithm. 
The  times  for  NDCOUNT  indicate  the  cost  of  visiting  every 
node  in  the  tree  exactly  once.  Hence  the  time  is  relatively 
constant  for  each  map  because  so  little  processing  is  done 
at  each  node.  This  value  also  gives  an  indication  of  the 
reliability  of  the  system  timing  routine  used.  Substan¬ 
tially  more  calculation  is  performed  by  AREA  with  more  vari¬ 
ation  with  respect  to  the  amount  of  time  spent  at  a  black 
node  vs.  a  gray  or  white  node.  The  conversion  frown  area  in 
pixels  to  area  in  acres  was  calculated  based  on  .142  acres 
per  pixel.  The  value  is  given  in  hundredths  of  an  acre, 
although  the  pixel  size  is  about  one  seventh  of  an  acre. 
The  coordinates  used  for  the  centroid  are  based  on  the 
upper  left-hand  corner  being  (0,0)  and  the  number  of  pixels 
in  both  directions  range  from  0  to  512.  The  same  coor¬ 
dinate  system  is  used  in  the  other  tables. 

The  REGION  PROPERTY  RESULTS  (Tables  5.5-7)  show  the 
cost  of  two  statistics  gathering  programs:  PERIMETER  and 
HANDW.  The  perimeter  is  measured  in  pixel  widths.  The 
enclosing  rectangle  calculated  by  HANDW  is  given  by  the 
coordinates  of  its  upper  left-hand  corner  and  its  width  and 
height.  HANDW  is  another  algorithm  that  treats  each  node 
equally  and  hence  produces  little  variation  in  its  timings 
within  a  given  map.  This  is  quite  different  from  PERIMETER, 
which  performs  four  FIND_NEIGHBOR  operations  for  each  black 
node;  hence  the  variations  in  the  cost  of  PERIMETER. 


The  data  for  the  WINDOW  program  is  presented  in  the 


WINDOW  RESULTS  (Tables  5.8-10).  The  window  used  is  the 
smallest  square  whose  width  is  a  power  of  two,  that  encloses 
the  smallest  bounding  rectangle  of  the  polygon  (calcu¬ 
lated  by  HANDW) ,  and  sharing  the  same  upper  left-hand 
corner  with  that  rectangle.  The  relation  between  the 
times  and  the  input  is  complicated,  as  it  is  effected  by 
both  the  size  of  the  window  and  the  greatest  common 
denominator  of  the  tree  size  and  the  two  coordinates  of  the 
upper  left-hand  corner.  The  smaller  the  greatest  common 
denominator  of  these  three  numbers,  the  greater  the  possible 
fracturing  of  large  nodes  in  the  input  tree. 

The  next  table,  INTERSECTION  STATISTICS  (Table  5.11), 
is  the  only  table  showing  a  binary  relation,  that  of  INTER¬ 
SECTION.  Three  large  regions  (the  center  of  the  flood 
plain  and  the  two  lowest  contours)  are  chosen  to  be 
intersected  with  the  land-use  classes  because  they  are  most 
likely  to  yield  interesting  results.  Since  the  center 
of  the  flood  plain  is  not  equivalent  to  a  contour  class, 
it  is  also  intersected  with  each  of  the  contour  classes. 
Note  that  the  cost  of  INTERSECTION  can  be  less  than  the  cost 
of  doing  a  NDCOUNT  on  both  trees  because  a  large  white 
node  in  one  tree  can  make  it  unnecessary  to  process  a  large 
subtree  in  the  other  tree.  As  well  as  the  cost  of  per¬ 
forming  the  INTERSECTION  operation  (measured  in  seconds) , 
the  table  also  gives  the  area  and  number  of  nodes  in  the 
result.  Note  that  a  UNION  table  is  not  shown  because  UNION 
behaves  in  the  same  manner  as  INTERSECTION  on  the  logical 
complement  of  the  inputs  (i.e.,  switch  the  black  and  white 
node  colors) .  Note  that  the  INTERSECTION  algorithm  is 
greatly  simplified  by  the  digitization  process's  alignment 
of  the  maps  so  that  the  pixel  at  (0,0)  corresponds  to  the 
same  ground  truth  in  each  map. 

The  final  table,  QUADTREE  TRUNCATION  STATISTICS  (Table 
5.12) ,  shows  the  amount  of  compression  one  can  obtain  by 
truncating  the  quadtree  maps.  The  usability  of  the  truncated 
quadtrees  is  discussed  at  the  end  of  Section  5.2  and  shown 
in  Figures  5. 4-5. 8.  Under  each  map's  name  there  are  two 
columns.  The  first  column  shows  the  number  of  nodes  in  the 
quadtree  that  is  formed  by  truncating  the  full  (depth  10) 
quadtree  to  the  depth  indicated  in  the  far  left  column.  The 
second  column  shows  the  percentage  of  nodes  in  the  full 
(depth  10)  quadtree  that  would  not  be  needed  for  the  trun¬ 
cated  quadtree. 
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TOPOGRAPHY  AREA  RESULTS 
(2  Or'  2) 
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(1  OF  5) 
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11 
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1 
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173 

3 

21 

1 .91 
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TABLE  5.7.  FLOODPLAIN  REGION  PROPERTY  RESULTS 


|  POLY-|PER-|TI*>jE|  ENCLOSING  |  TIMS  I 
I  GON  |  Iti-|  IN  |  RECTANGLE  j  IN  | 
I  | ETEK | SECS i  X  Y  WID  hGT j SECS | 


I  right  1 1776 1  0.91105  0  295  450 1  0.4 1 

I  left  |2270|  0.91  0  0  274  450|  0.4| 

I  center  |  2642  |  1  .0 1  3  0  280  450 1  0.4| 
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TABLE  5 .fa  .  LAN DUSE  WINDOW  RESULTS 


U 

of  5) 

I  POLY- I 

1  GON  | 

WINDOW 

EX  |  FY  | WIDTH 

1  TIME  1 
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LAN DUSE  WINDOW  KE6ULT8 


(2 

of  5) 
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1  G ON  | 

WINDOW 

FX  |  FY  I  WIDTH 
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LAW  DU  Si.  HlWDOU  RLSULT6 
(3  or  5) 
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1  3.9 

l  r  .3  1 

0| 

0| 

512 

1  39  .4 

1  r.4  1 

2601 

ol 

64 

1  1.5 

1  r.5  I 

2731 

46 1 

126 

111  .5 

1  ucb  .1  1 

70 1 

41 

32 

1  0.4 

lucb.2  1 

44  j 

99 1 

16 

1  0.3 

lucc.l  1 

46 1 

4 1 

32 

1  0.3 

luce. 2  1 

511 

501 

16 

1  0.3 

luce. 3  | 

391 

74| 

16 

1  0.3 

luce  .4  1 

74 1 

701 

6 

1  0.0 

luce. 5  | 

621 

771 

1© 

I  0.4 

1 ucc  .6  1 

1751 

3731 

32 

I  1  .4 

LANDUSE  WINDOW  RESULTS 
(4  of  5} 
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i 

7.4  1 

302| 

1421 

4 1  0.01 

i 

7.5  | 

3641 

1401 

1)  O.OI 

i 

6.1  i 

1631 

301 

64)  5.7) 

i 

8.2  1 

3101 

581 

8 1  O.li 

i 

6.3  | 

2561 

53  1 

256139.6) 

i 

6.4  1 

3401 

51 1 

2 1  0.0) 

i 

8.5  1 

3651 

66 1 

1261  4.6| 

i 

6.6  j 

3301 

1911 

2 1  0.0) 

i 

9.1  | 

1921 

341 

641  3.2| 

i 

9.2  1 

309] 

991 

16)  0.4) 

i 

9.3  i 

3261 

67) 

321  1.1| 

i 

9.4  | 

325| 

701 

126111 .31 

i 

9.5  1 

3611 

1951 

4 1  O.li 

i 

9.6  | 

3721 

2021 

161  O.li 

i 

10.1  i 

195 1 

40 1 

321  1.5) 

i 

10.2  1 

3671 

76  l 

161  0.4( 

i 

10.3  | 

3501 

951 

6 1  0.1) 

i 

10.4  | 

3341 

1061 

1261  5.11 

i 

11.1  | 

347| 

1131 

16  1  0  .4  | 

i 

11.2  1 

362) 

173) 

4 i  0.0) 

TABLE  5.10.  FLOODPLAIN  WINDOW  RESULTS 

I  POLY— i  WINDOW  I TIME  I 

|  GUN  |  FX  |  FY  | WIDTH | SECS  I 


149 


I 


TABLE  5.11.  INTERSECTION  ST  AT  IS  VC  i> 

(1  Of  3) 


i  1 

I  TREE  1  I 

1 

TREE  2  1 

AREA  ! 
PlXELSl 

AREA  | 
ACRES  | 

NUMBER  OF  NODES  | 
GRAY  1  BLACK | WHITE  1 

TXNEl 
SECS  I 

1 f  .center  I 

t.l  1 

28446 ! 4039  .33  I 

1672| 

23941 

26231 

5.0 1 

1 f  .center  1 

t  .2  1 

1 2fal  | 

181  .901 

7981 

7601 

16151 

2.6| 

1  f  .center  i 

t.3  j 

ol 

o.ool 

ol 

0| 

1 1 

0.1 1 

1  f  .center  1 

t.4  1 

ol 

0  .001 

Ol 

ol 

1 1 

o.ol 

1  f  .center  i 

t.5  | 

ol 

o.ool 

Ol 

0| 

1 1 

o.ol 

i f  .center  1 

t  .6  | 

ol 

0.001 

Ol 

0| 

1 1 

o.ol 

1  f  .center  1 

t  .7  | 

ol 

o.ool 

Ol 

ol 

11 

o.ol 

I  f  .center  1 

t  .8  | 

ol 

o.ool 

ol 

0| 

1| 

o.ol 

1  £  .center  1 

t  .9  | 

ol 

o.ool 

ol 

ol 

1| 

o.ol 

i  t  .center  i 

t  .10  1 

ol 

o.ool 

ol 

0| 

li 

o.ol 

1 f  .center i 

t  .11  j 

ol 

o.ool 

Ol 

0| 

1 1 

o.ol 

1  f  .center  1 

l.acc  1 

14721 

209  .021 

2231 

2761 

3921 

0.7| 

1  f  .center  1 

l.acp  i 

1521 

21 .581 

721 

66| 

1491 

0.31 

1  f  .center  1 

1  .ar  j 

01 

o.ool 

Ol 

0| 

1 1 

o.ol 

1 1  .center  i 

1 .are  i 

ol 

o.ool 

ol 

0| 

1 1 

o.ol 

1  f  .center  1 

1  .avf  | 

58691 

833 .401 

12251 

15551 

2121 1 

4.0| 

1  f  .center  1 

1  .aw  1 

11376|1615.39| 

12771 

16531 

21791 

4.11 

1  f  .center  1 

1  .bbr  1 

4321 

61  .34| 

134| 

1531 

2501 

0  .4  | 

i  £  .center  1 

l.beq  1 

229  1 

32 .521 

88) 

97 1 

1681 

0.31 

1  f  .center  I 

1  .be  s  1 

1471 

20.871 

48 1 

SI  | 

94 1 

0.1 1 

1  f  .center  1 

l.bt  | 

1321 

18.741 

46 1 

511 

601 

0.1 1 

I  f  .center  1 

l.fo  1 

4691 

66.601 

1661 

170| 

321| 

0.5| 

1  f  .center  1 

l.lr  | 

9051 

128.511 

3541 

4161 

6471 

1  .ol 

1  f  .center  1 

l.r  1 

46 1 

6.531 

421 

341 

931 

0.1  I 

1  f  .center  1 

1 . ucb  | 

ol 

o.ool 

ol 

ol 

11 

o.ol 

1  f  .center  1 

1 . ucc  1 

3  | 

0.431 

131 

3| 

371 

0.1 1 

1  f  .center  1 

1 . ucr  j 

ol 

o.ool 

ol 

0| 

ll 

o.ol 

1  f  .center  1 

1 .  ucw  | 

ol 

o.ool 

Ol 

0| 

1 1 

o.ol 

1  f  .center  1 

1  .vies  | 

12861 

182.611 

2751 

3351 

491 1 

1  .01 

1 f  .center  j 

l.uil  | 

ol 

o.ool 

Ol 

o| 

1 1 

o.ol 

1  £  .center  1 

1  .uis  | 

107! 

10.651 

35 1 

381 

66 1 

0.1  I 

1  £  .center  1 

1 .  uiw  | 

Ol 

0  .001 

Ol 

ol 

1 1 

o.ol 

1  £  .center  i 

1 .  unk  | 

ol 

o.ool 

ol 

ol 

1 1 

o.ol 

1  £  .center  I 

1 . uoc  1 

ol 

o.ool 

0| 

0| 

ll 

o.ol 

1  f  .center  j 

1  .uog  1 

oj 

o.ool 

ol 

ol 

11 

o.ol 

1  £  .center  1 

1  • uoo  | 

751 

10.651 

31 1 

331 

61 1 

0.1 1 

I  £  .center  1 

1 .  uop  i 

1841 

28.131 

471 

55| 

871 

0.21 

1  £  .center  1 

1 . uov  | 

ol 

o.ool 

ol 

0| 

1 1 

o.ol 

!  £  .center  1 

1 .  urh  | 

201 

2  .841 

19  1 

111 

471 

0.1 1 

1  £  .center  j 

1  .urs  | 

21801 

309  .561 

6171 

9351 

15171 

2.51 

i  £  .center  I 

1 . uus  i 

2491 

35.361 

651 

72| 

124| 

0.21 

i  £  .center  1 

1 . uut  1 

2241 

31 .81 1 

77| 

83| 

1491 

0.31 

i  £  .center  1 

1  .w  | 

108| 

15 .341 

38 1 

39  | 

76 1 

0  .2  i 

1  £  .center  I 

1  .wo  1 

661 1 

93 .86| 

1211 

1391 

2251 

0.31 

1  £  .center  1 

1  »W8  | 

3401! 

482 .941 

11701 

14611 

20301 

3.4| 

1  £  .center  1 

1  .wwp  j 

ol 

o.ool 

Ol 

ol 

ll 

o.ol 
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t  .1 
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t  .1 
t  .1 
t  .1 
t  .1 
t  .1 
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t  .1 
t  .1 
t  .1 
t  .1 
t  .1 
t.l 
t  .1 
t.l 
t.l 
t.l 
t.l 
t.l 
t.l 
t.l 
t  .1 
t.l 
t.l 
t  .1 


INTERSECTION  STATISTICS 
(2  of  3) 


1 

1 

TKEE  2 

I  AREA  | 

1 PIXELSl 

AREA  | 
ACRES  | 

NUMBER  OF  NODES  | 
GRAY | BLACK | WHITEl 

TIME 

SECS 

1 

1  .acc 

1 

39071 

554.791 

6561 

844| 

11251 

2.0 

1 

1  .acp 

1 

1402  i 

199 .081 

408 1 

4901 

7351 

1  .2 

1 

1  .ar 

1 

581  | 

82.501 

1641 

2031 

2901 

0.5 

1 

1  .are 

1 

oi 

o.ool 

ol 

ol 

1 1 

0.0 

1 

1  .avf 

1 

163581 

2322 .841 

21071 

29721 

33501 

6.5 

1 

1  .aw 

I 

202881 

2880.901 

2180| 

29841 

3557! 

6  .6 

1 

1  .bbr 

1 

4321 

61 .34| 

134| 

1531 

250| 

0.4 

1 

1  .beq 

1 

229  | 

32 .521 

88 1 

97| 

168| 

0.3 

1 

1  .bes 

1 

114| 

15  .191 

52  | 

54 1 

1031 

0.2 

1 

l.bt 

1 

88 1 

12  .501 

34 1 

37  | 

66 1 

0.1 

1 

l.fo 

1 

3811 

54  .101 

1651 

1681 

328| 

0.5 

1 

l.lr 

1 

9131 

129  .651 

3591 

424| 

6541 

1  .1 

1 

1  .r 

1 

25| 

3.55| 

381 

25| 

90 1 

0.1 

1 

1 .  ucb 

1 

01 

o.ool 

ol 

ol 

1 1 

0.0 

1 

1 .  ucc 

1 

139  i 

19  .74| 

531 

58 1 

1021 

0.2 

1 

1 .  ucr 

1 

7691 

109  .20| 

1791 

2231 

3151 

0.5 

1 

1 .  ucw 

1 

3051 

42.311 

1121 

140| 

197  i 

0.3 

I 

1 .  ues 

1 

14041 

199  .371 

348| 

4471 

5981 

1  .0 

1 

1  .uil 

1 

371| 

52.681 

89  | 

101 1 

1671 

0.3 

I 

1  .uis 

1 

6271 

89.031 

1771 

204| 

328| 

0,5 

1 

1  .uiw 

1 

oi 

o.ool 

ol 

ol 

11 

0.0 

1 

1  .unk 

1 

01 

o.ool 

ol 

ol 

1| 

0.0 

1 

1  .uoc 

1 

ol 

o.ool 

0! 

ol 

li 

0.0 

1 

1 .  uog 

1 

01 

o.ool 

ol 

ol 

l! 

0.0 

1 

1 .  uoo 

1 

4901 

69  .581 

1071 

121| 

201 1 

0.3 

1 

1  .uop 

1 

148 1 

21 .021 

50 1 

58 1 

931 

0.2 

l 

1  .uov 

1 

ol 

o.ool 

0| 

ol 

1 1 

0.0 

1 

1 .  urh 

1 

1261 

17.891 

321 

331 

64 1 

0.1 

1 

1 .  urs 

1 

38511 

546.841 

1336| 

15081' 

25011 

4.0 

1 

1 .  uus 

1 

2611 

3? .06| 

74 1 

81 1 

1421 

0.2 

I 

1 .  uut 

1 

938| 

133  .201 

314| 

3591 

5841 

1  .0 

1 

1  .vv 

1 

1 08  | 

15  .341 

38 1 

391 

76 1 

0.1 

1 

1  .wo 

1 

6611 

93  .86| 

1211 

1391 

2251 

0.4 

1 

1  .ws 

1 

34021 

483 .081 

11681 

14821 

20231 

3.4 

1 

1  .wwp 

1 

ol 

o.ool 

ol 

ol 

1 1 

0.0 
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INTERSECTION  STATISTICS 
(3  of  3) 


i 

1 

TREE  1 

I  AREA  i 
TREE  2  | PIXELSl 

AREA  l 
ACRES  | 

NUMBER  OF  NODES  I 
GRAY | BLACK ! WHITE 1 

TIME  I 
SECS  | 

1 

t  .2 

1  .acc  I 

2434! 

345  .63! 

506! 

616 1 

9031 

1  .51 

I 

t  .2 

l.acp  j 

12652 ! 

1624 .96 | 

17501 

2412! 

26391 

5.5! 

1 

t  .2 

1  .ar  i 

3401 

46.261 

130i 

1301 

261  i 

0  .4 1 

1 

t  .2 

l.are  j 

1521 

21 .38! 

32! 

261 

71| 

0.1 1 

1 

t  .2 

l.avf  i 

6716! 

953.671 

1163! 

1469| 

20211 

3.81 

i 

t  .2 

1  .aw  j 

7922  I 

1124.92! 

15351 

1886! 

27201 

4.9  1 

I 

t  .2 

l.bbr  1 

oi 

o.ooi 

o! 

ol 

1 1 

0.0! 

1 

t  .2 

1  .beq  j 

Ol 

0  .00! 

0! 

ot 

1 1 

o.oi 

1 

t  .2 

1  .be s  j 

33! 

4.69| 

25 1 

211 

551 

0.1 1 

1 

t  .2 

l.bt  | 

941 ! 

133 .621 

2261 

257| 

428! 

0.7! 

1 

t  .2 

l.fo  1 

10091 

143 .261 

3451 

3941 

6421 

1  .0! 

1 

t  .2 

1  .lr  | 

35| 

4.97| 

42! 

29| 

98 1 

0.21 

1 

t  .2 

l.r  I 

2209! 

313  .67! 

443 1 

5711 

759| 

1  .4! 

i 

t  .2 

1 . ucb  ! 

249  i 

35  .36! 

62! 

69| 

1161 

0.2| 

1 

t  .2 

1 . ucc  | 

679  | 

124.821 

1721 

1831 

334| 

0.5! 

1 

t  .2 

1 . ucr  | 

749! 

106.361 

1791 

2241 

314| 

0.6! 

1 

t  .2 

1 .  ucw  | 

ol 

o.ooi 

o! 

ol 

1| 

o.ol 

1 

t  .2 

1 . ues  | 

224! 

31  .81 1 

106! 

1131 

206| 

0  .3  | 

1 

t  .2 

l.uil  | 

51! 

7.241 

37| 

33! 

791 

0.1  i 

1 

t  .2 

l.uis  j 

415! 

58.931 

1371 

1511 

261! 

0.4| 

1 

t  .2 

1 .  uiw  j 

161  | 

22.861 

69! 

71! 

1371 

0.21 

1 

t  .2 

1  .unk.  | 

6601 

93.721 

99  i 

1261 

1721 

0.3| 

1 

t  .2 

1 . uoc  1 

262! 

40.041 

46  { 

571 

821 

0.1  i 

1 

t  .2 

1 . uog  | 

7651 

108.63! 

1401 

1951 

226| 

0  .4 1 

1 

t  .2 

1 . uoo  | 

ol 

0.00! 

ol 

ol 

1 1 

o.oi 

1 

t  .2 

1 . uop  j 

65! 

9.231 

51 ! 

321 

1221 

0.2! 

1 

t  .2 

1 . uov  j 

171| 

24.26! 

57| 

721 

1001 

0  .2  | 

1 

.  t.2 

1  .urh  | 

41 1 

5.821 

28| 

201 

651 

0.1 1 

1 

t  .2 

1 .urs  j 

15611 1 

2245  .161 

20791 

28661 

33721 

6  .41 

I 

t  .2 

l.uus  | 

ol 

o.oot 

0| 

o! 

1 1 

0.01 

1 

t.2 

1 . uut  I 

961  t 

139 .301 

491 1 

570! 

9041 

1  .5! 

1 

t.2 

1  .vv  | 

0| 

o.oo! 

ol 

o! 

1 1 

o.oi 

1 

t  .2 

1  .wo  j 

ol 

0.00! 

o! 

0) 

11 

o.oi 

i 

t  .2 

1  .ws  | 

7  | 

0.99| 

29  | 

71 

611 

0.3! 

1 

t  •  2 

1  .wwp  j 

116! 

16  .761 

63! 

551 

1351 

0.2! 

TABLE  5.12.  QUADTREE  TRUNCATION  STATISTICS  FOR  EACH  MAP 


1 

DEPTH 

I  LANDUSE  MAP 

1 

TOPOGRAPHY  MAP | 

FLOODPLAIN  MAPI 

1 

OF 

|  NUM  OF | 

%  RE- 

1 

NUM  OF | 

%  RE- 

1 

NUM  OF | 

%  RE- 

1 

1 

TREE 

|  NODES 

i 
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1 
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1 

DUCED 

1 
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1 
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1 

1 

10 

i  3B233 

1 

00.00 

1 

33349 

1 
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\ 

6941 

1 

00.00 

l 

1 

9 

1  220B9 

1 

44.22 

1 

18517 

1 

44.47 

I 

4473 

1 

35.55 

1 

1 

b 

1  9489 

1 

75  .18 

I 

7473 

1 

77.59 

1 

2297 

1 

66.91 

1 

1 

7 

I  3341 

1 

91  .26 

1 

2537 

1 

92.39 

1 

1093 

1 

84.26 

1 

1 

6 

1  1057 

1 

97  .23 

1 
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1 

97.50 

1 
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1 

92.38 

1 

1 

5 

I  309 

1 

99  .19 

I 

296 

1 
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1 

213 

1 

96  .94 

1 

1 

4 

1  85 

1 

99  .78 

1 

77 

1 

99  .77 

1 

77 

1 

98.89 

1 

1 

3 

I  21 

1 

99  .95 

1 

21 

1 

99  .94 

1 

21 

1 
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1 

J 

2 
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1 

99  .99 

I 

5 

1 

99  .99 

5 

1 

99  .93 

1 

1 

1 

1  1 

1 

99  .99 

1 

1 

1 

99  .99 

i 

1 

1 

99  .99 

1 
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7.  Conclusions  and  future  plans 
7.1.  Conclusions 


This  project  gave  a  firm  empirical  basis  to  much  of  the 
theoretical  analysis  previously  undertaken  for  quadtrees 
both  as  to  their  structure  and  their  algorithmic  efficien¬ 
cies.  In  particular,  the  following  conclusions  should  be 
noted: 

(1)  Errors  in  the  calculations  of  properties  encoded  by 
quadtrees  (e.g.,  areas  and  perimeters  of  various  land 
use  classes)  are  due  entirely  to  errors  introduced  by 
the  original  digitization.  No  new  errors  are  introduced 
by  quadtree  manipulation. 

(2)  Significant  reductions  in  file  size  are  achieved  when 
an  image  is  converted  from  a  binary  array  representa¬ 
tion  to  a  quadtree  representation.  This  is  true  for 
both  the  multicolored  and  black/white  cases. 

(3)  The  block  decomposition  of  the  image  resulting  from  the 
quadtree  representations  yields  major  increases  in 
display  speed. 

(4)  Truncation  of  quadtrees  can  be  used  to  generate  reason¬ 
able  image  approximations  that  are  consistently  more 
compact . 

(5)  Quadtree  algorithms  are  easy  to  implement  in  structured 
programming  languages  (e.g.,  C) . 

(6)  Neighbor  finding  was  found  to  require  visiting  3.5 
nodes  on  the  average  for  each  instantiation.  This  was 
even  better  than  what  was  expected  theoretically. 

(7)  Ropes  (an  alternative  neighbor  finding  technique)  were 
found  to  be  not  worth  the  added  expense  of  extra 
storage. 

(8)  Set  operations  such  as  union  and  intersection  are  effi¬ 
cient  and  can  be  used  to  extract  information  from 
images  containing  different  properties. 

It  should  also  be  noted  (in  conjunction  with  (3)  and  (4) 
above)  that  quadtrees  could  be  used  effectively  in  image 
transmission,  enabling  the  viewer  to  recieve  a  very  compact 
approximation  of  the  image  followed  by  a  series  of  modifica¬ 
tions  that  render  the  image  increasingly  more  precise. 

7.2.  Future  plans 


The  first  phase  of  this  project  has  dealt  with  digiti¬ 
zation  of  a  government-furnished  geographic  database  and  its 


representation  in  quadtree  form;  and  with  development  of 
algorithms  for  basic  operations  on  quadtree-represented 
regions  (set-theoretic  operations,  point- in- region  determi¬ 
nation,  region  property  computation;  submap  generation) .  The 
efficiency  of  these  algorithms  was  studied  theoretically  and 
experimentally. 

The  following  tasks  are  planned  for  the  second  phase: 

(a)  Query  language.  Design  of  a  high-level  query  language 
permitting  easy  interaction  with  the  database  by  users, 
thus  making  the  quadtree  representation  transparent  to 
the  users. 

(b)  Database  updating .  Develpment  of  algorithms  for  addi¬ 
tion^  deletion,  and  editing  of  data  items  in  a 
quadtree-encoded  database. 

(c)  Point  and  linear  feature  data.  Quadtree-like  data 
structures  will  also  Ee  used  for  the  storage, 
retrieval,  and  editing  of  point  geographic  data.  Algor¬ 
ithms  will  be  incorporated  for  performing  these  func¬ 
tions  and  for  interfacing  between  tree  representations 
of  point  and  area  data.  Recently,  quadtree-like  data 
structures  have  been  developed  for  representing  region 
borders  and  curves.  The  interface  between  these  struc¬ 
tures  and  the  tree  representations  of  points  and 
regions  will  be  investigated. 
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Appendix:  Facilities  used 


Two  computers  produced  by  the  Digital  Equipment  Cor¬ 
poration  are  used  by  this  project.  Program  development  and 
small-scale  testing  are  performed  on  a  PDP  11/45.  Our  PDP 
11/45  has  a  256k  bytes  of  actual  memory  of  which  only  64k 
bytes  are  directly  addressable,  no  virtual  memory  capabili¬ 
ties,  a  disk  fetch  speed  of  1.2  megabits/second,  and  a 
memory  cycle  speed  of  approximately  500  microseconds.  The 
execution  times  in  the  tables  of  this  report  refer  to  the 
execution  speed  on  the  VAX  11/780.  The  VAX  11/780  has  2000k 
bytes  of  actual  memory,  6000k  bytes  of  virtual  memory,  a 
disk  fetch  speed  of  approximately  0.6  megabits/second,  and 
a  memory  cycle  speed  of  approximately  1400  nanoseconds. 
The  size  of  a  quadtree  node  is  12  bytes  on  the  PDP  11/45 
and  24  bytes  on  the  VAX  11/780.  This  difference  is  caused 
by  the  different  word  size  on  each  machine.  Both  the  PDP 
11/45  and  the  VAX  11/780  run  the  UNIX  operation  system 
(versions  6  and  7  respectively) . 

The  picture  output  device  used  by  this  project  is  a 
Grinnell  GMR-27  Display  Processor.  Its  memory  consists  of 
thirteen  512x512  bitplanes.  Twelve  of  these  bitplanes  carry 
color  information  (4  bits  for  each  of  the  colors:  blue, 
green,  and  red) .  The  thirteenth  bitplane  is  used  for  a 
white  overlay  capability.  The  high  order  eight  bitplanes  of 
the  twelve  color  bitplanes  can  also  be  displayed  to  create  a 
grayscale  output.  The  output  speed  of  quadtrees  on  this 
device  is  considerably  faster  than  a  raster  scan  output  of  a 
picture  file,  because  the  GMR-27  can  output  a  rectangle  on 
the  display  screen  directly  from  the  rectangle's  coordinates 
(i.e.,  a  separate  command  is  not  necessary  for  each  pixel  in 
the  rectangle  as  is  done  when  a  picture  file  is  output  in 
raster  scan  mode) . 


As  our  display  device  is  connected  to  a  computer  with 
restricted  memory  (see  Appendix),  we  will,  in  addition  to 
the  above,  be  investigating  more  compact  in-core  representa¬ 
tions  and  the  effect  of  user-controlled  paging  on  algorithm 
efficiencies.  This  will  be  done  in  conjunction  with  the 
development  of  a  quadtree  editor  (which  requires  interactive 
use  of  display  device) . 


